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ABSTRACT 

Mathematical  models  are  developed  for  the  hydrodynamic 
and  aerodynamic  forces  which  act  on  a  sailing  boat.   The 
solution  of  the  force  and  moment  equilibrium  provides  the 
sailing  boat's  speed  and  attitude  for  performance  prediction. 
Each  force  model  is  evaluated  for  its  capability  to  make 
adequate  engineering  predictions.   Further,  expressions  are 
given  for  the  dimensions  and  coefficients  to  provide  a  feasible 
baseline  design  to  which  parametric  variations  may  be  made. 

The  system  of  force  equations  constitute  a  parameterized 
system  model  in  four  degrees  of  freedom  for  steady  state 
sailing  conditions.   The  mathematical  model  provides  a  design 
tool  for  the  naval  architect  conducting  geometric  parameter 
trade-off  analysis  against  hydrodynamic  performance. 

The  principal  result  is  a  computerized  model  which  is 
capable  of  predicting  boat  speed,  heel  angle,  leeway  angle 
and  rudder  angle  for  a  geometrically  defined  sailing  boat  in 
a  given  wind  condition.  'Example  boats  are  evaluated  and  the 
results  recorded  in  order  to  demonstrate  the  use  and  versa- 
tility of  the  model. 
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NOMENCLATURE 


A,  =  keel  cross-sectional  area. 
k 

A,  =  projected  lateral  area  of  the  hull  above  the  LWL. 

A  =  cross-sectional  area  of  the  maximum  hull  station, 


m 


or  cross-sectional  area  of  mast  depending  on 
context. 


A  =  platform  area  of  the  rudder, 

r  ^ 

A  =  total  area  of  the  sails. 
s 

A  =  projected  transverse  area  of  the  hull  above  the  LWL 

BWL  =  maximum  waterline  beam. 

C-,  =  rudder  cross  flow  drag  coefficient. 

C,  =  rudder  section  form  drag  coefficient, 
do                              3 

C,.  =  sail  induced  drag  coefficient. 

Cp  =  friction  drag  coefficient. 

Cpfh  =  fairbody  friction  drag  coefficient. 

Cp,  =  keel  friction  drag  coefficient. 

Cf  =  hull  form  drag  coefficient. 

Cj.  =  rudder  friction  drag  coefficient. 

fr  3 

C,  =  rudder  lift  coefficient, 
lr 

C,  =  sail  lift  coefficient. 
Is 

C   -,.  =  rudder  moment  coefficient  about  the  mean  quarter 
mc/4        ,    , 

'  chord. 

C  =  rudder  normal  force  coefficient. 
n 

C  =  prismatic  coefficient. 

C  =  parent  prismatic  coefficient. 

C  =  optimum  prismatic  coefficient. 

C  =  rudder  drag  coefficient. 
xr                   3 
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C  =  keel  side  force  coefficient 

y 

C  i_  —  slender-body  side  force  coefficient, 
ysb  2 

C  =  rudder  side  force  coefficient. 
yr 

D  =  sail  drag 

D  =  off wind  sail  drag. 
so  r 

F,  =  aerodynamic  hull  force. 

F  =  Froude  Number. 
n 

H  =  maximum  canoe  body  draft. 

K  =  heeling  moment. 

K,  =  hydrostatic  righting  moment 

K,  =  keel  heeling  moment. 

K  =  rudder  heeling  moment. 

K  =  sail  heeling  moment. 

L  =  load  waterline  length. 

L  =  sail  lift. 
s 

L  =  offwind  sail  lift. 
so 

LCB  =  position  of  center  of  buoyancy  (%LWL) . 

LCB  =  parent  position  of  center  of  buoyancy. 

LCB  =  optimum  location  of  center  of  buoyancy. 

LWL  =  load  waterline  length. 

M  =  pitching  moment 

N  =  yawing  moment 

N,  =  keel  yawing  moment. 

N  =  rudder  yawing  moment. 

N  b  =  slender-body  yawing  moment. 

N  =  yawing  moment  due  to  sail  driving  force 

N  =  yawing  moment  due  to  sail  side  force. 
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P  =  vector  of  assumed  variable  values. 

p,  ~  -,  =assumed  value  of  individual  variable. 
1,2,3,  or  4 

Rf  =  total  frictional  resistance. 

R_  =  total  form  resistance, 
fm 

R  =  Reynold's  Number. 

R  =  upright  resistance. 

R  =  wave  resistance. 
w 

R  =  parent  wave  resistance. 

wp  e 

R  =  vector  of  force  equilibrium  errors. 

R,  0  ~.  .=   individual  force  equilibrium  errors. 

Sf,  =  fairbody  wetted  surface. 

S.  =  keel  wetted  surface. 
k 

S  =  rudder  wetted  surface, 
r 

T  =  canoe  body  maximum  draft. 

U,  =  boat  velocity  in  the  X  direction, 
b                     J                      o 

V  =  apparent  wind  velocity, 
a 

V,  =  boat  velocity  parellel  to  the  y  axis. 

V  =  true  wind  velocity. 

X  =  force  in  the  X   direction. 

o 

X,  =  hull  resistance. 
h 

X,  =  keel  and  hull  induced  drag. 

X  =  axis  on  the  plane  of  the  free  surface  in  the 
o 

direction  of  boat  velocity. 

X  =  rudder  induced  drag. 
r 

X  =  sail  driving  force. 

V  =  force  in  the  Y   direction,  side  force. 

Y,  =  keel  side  force. 
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Y  =  an  axis  on  the  plane  of  the  free  surface  at 

right  angles  to  boat  velocity. 

Y  =  rudder  side  force. 

r 

Y  =  sail  side  force. 
s 

Z  =  force  in  the  Z   direction,  heave  force. 

Z  =  an  axis  perpendicular  to  the  free  surface 
downward . 

af  =    .  fortriangle  aspect  ratio. 

a  =  mainsail  aspect  ratio, 
m 

a  =  rudder  aspect  ratio. 

a„  =  total  sail  area  aspect  ratio, 
s 

b,  =  draft  at  forward  perpendicular. 
b  =  maximum  draft. 

b  =  rudder  span. 

b-,  2  3=s  intermediate  draft  measurements. 

c  =  rudder  root  chord. 

c  =  rudder  tip  chord. 

rt  r 

c  ..  =  rudder  quarter  chord. 

c,  =  keel  mean  chord, 
k 

c  =  rudder  mean  chord. 

d1  =  metacentric  heigh t/LWL. 

d2  =  non-linear  hydrostatic  righting  moment  coefficient. 

d,  =  Froude  Number  dependent  righting  moment  coefficient , 
f .  (  )  =  spline  cubic  function. 

g  =  acceleration  due  to  gravity. 

r  =  canoe  body  draft  at  maximum  draft  measurement. 

s  =  sinkage. 


s.  •    =      9R./3P  •  • 
ID  i    D 
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t     =     taper  ratio. 

x     =     axis  in  the  load  waterline  plane  along  the 
centerline  of  the  boat. 

x,     =     horizontal  distance  to  the  maximum  draft 
b  measurement. 

x  .   =     horizontal  distance  to  the  aerodynamic  center 


'cph 


cpr 


z 


zcpxs 


of  effort  on  the  hull. 


x  ,   =     horizontal  distance  to  center  of  effort  on 
cpK        the  keel. 

x     =     horizontal  distance  to  the  center  of  effort  on 


the  rudder. 


x  =     horizontal  distance  to  center  of  effort  on 
cp         the  sail. 

x  =     main  boom  length. 

x.  =     base  of  the  fortriangle. 
3 

x  =     horizontal  distance  to  the  mast, 
m 

x  =     0.0,  location  of  maximum  draft. 

x  -  /A  -      horizontal  distance  from  the  rudder  mean  chord 
*  quarter  point  to  the  rudder  center  of  effort. 

x  =     horizontal  distance  to  the  rudder  stock, 
rs 

X,  2  -=     horizontal  distance  from  maximum  draft  to 
'  '       intermediate  draft  measurements. 

y  =     axis  orthaginal  to  the  x  axis  to  starboard 
in  the  waterline  plane. 

z  =     axis  orthaginal  to  the  waterline  plane,  downward, 

z   .  =     vertical  distance  to  the  center  of  effort  on 
p         the  keel. 


=     vertical  distance  to  the  center  of  effort  on 
p         the  rudder. 

vertical  distance  to  the  sail  driving  force. 

vertical  distance  to  the  sail  side  force. 


zcpys 

z£u   =     freeboard. 
fb 
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z  =  vertical  distance  to  the  center  of  gravity. 

g 

z.  =  mast  height. 

z  =  main  hoist. 
P 

z  =  vertical  distance  from  the  rudder  root  chord  to 


ro 


the  center  of  effort  on  the  rudder. 


z  =  vertical  distance  to  the  rudder  root  chord, 
rrc 

V  =  volumetric  displacement  of  the  fairbody. 
V,  =  volumetric  displacement  of  the  keel. 

A  =  weight  of  the  boat  in  tons. 

a  =  rudder  angle. 

8  =  apparent  wind  angle. 

8  =   _   upright  apparent  wind  angle. 

Y  =  true  wind  angle. 

6  =  variable  increment. 

0  =  pitch  angle. 

A  =  angle  between  boat  centerline  and  aerodynamic 
hull  force. 

Ytt  r.  =  kinematic  viscosity  of  v/ater . 

H2° 

it  =  3.14159 

p  .  =  density  of  air. 
air 

Ptr  r>  =  density  of  water. 

H2° 

<J)  =  heel  angle. 

ip  =  leeway  angle. 

ft  =  rudder  sweep  angle. 
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I.   INTRODUCTION 


1.1  General 


All  designers  of  engineering  systems  desire  a  procedure 
for  predicting  the  performance  of  their  designs  prior  to 
final  construction.   Until  lately,  prediction  of  a  sailing 
yacht's  performance  has  been  limited  to  evaluation  of  scale 
models  and  intuition.   The  development  of  the  high  speed 
computer  and  numerical  techniques  for  rapid  solution  of  non- 
linear equations  has  made  the  use  of  mathematical  models  a 
reasonable  alternative  for  performance  prediction.   It  is 
the  purpose  of  this  paper  to  present  mathematical  models  for 
the  forces  and  moments  acting  on  a  sailing  boat.   The  solution 
of  the  force  and  moment  equilibrium  will  provide  the  boat's 
speed  and  attitude  for  performance  prediction. 

Each  force  is  parameterized  in  dimensions  and  coefficients, 
generally  utilized  by  designers,  to  make  the  model  a  useful 
design  tool  for  the  naval  architect.   Further,  expressions 
are  given  for  the  required  dimensions  and  coefficients  to 
provide  a  feasible  baseline  design  from  which  parametric 
variations  may  be  made.   Each  force  model  is  evaluated  for 
its  capability  to  make  adequate  engineering  predictions, 
including  inherent  limitations. 

1.2  Historical  Development 

From  the  earliest  days,  yacht  design  has  been  called  an 

"art".   Though  the  use  of  wind  propulsion  for  ships  is  ancient, 

the  action  of  the  forces  and  moments  on  the  sailing  ship  is 
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complex  and  multidimensional.   Historically  there  has  been 
no  technique  for  simple  evaluation  of  the  effect  a  change  in 
one  parameter  has  on  the  forces  acting  on  the  boat. 

Until  fifty  years  ago  design  evolved  slowly  through  ober- 
vation  of  existing  ships.   Due  to  the  inability  to  evaluate 
a  change  prior  to  incorporation  into  a  newly  constructed 
yacht,  new  designs  tended  to  consist  of  only  conservative 
changes  to  a  previously  successful  design.   This  was  under- 
standable for  failure  of  an  innovation  in  the  final  yacht 
could  be  professionally  catastrophic  for  the  designer. 

In  the  early  1900' s,  techniques  were  developed  for  pre- 
dicting yacht  performance  from  model  tests.   A  designer  could 
now  build  a  relatively  inexpensive  model  and  evaluate  it 
against  other  model  results  prior  to  incorporating  a  modifica- 
tion into  a  new  design.   Adequate  prediction  of  full  scale 
results  have  been  limited  by  adverse  scale  effects,  but  the 
usefulness  of  these  tests  is  unquestioned.   Unfortunately, 
the  cost  of  model  testing  still  precludes  extensive  parametric 
trade-offs.   Each  parametric  change  requires  construction  and 
testing  of  individual  models.   Usually,  the  limited  design 
budget  for  one  yacht  precludes  any  systematic  variation  of 
parameters,  and  testing  is  limited  to  one  or  two  candidate 
models.   The  appeal  of  a  parameterized  mathematical  model 
which  is  inexpensive  to  utilize  is  apparent. 
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The  idea  of  using  mathematical  expressions  for  describing 
yacht  performance  is  by  no  means  a  new  concept.   The  rating 
rules  used  to  handicap  yachts  are  an  attempt  to  apply 
mathematical  models  to  boat  performance.   Because  rating 
rules  attempt  to  measure  a  boat's  speed  potential  for  a  large 
range  of  conditions  and  reduce  the  results  to  a  single  number, 
it  is  easily  understood  why  they  have  little  application  to 
predicting  performance  for  a  given  wind  direction  and  speed. 

It  was  not  until  the  second  half  of  this  century  that 
towing  tank  test  results  were  combined  with  sail  force 
coefficients  in  a  computer  to  solve  force  and  moment  equili- 
briums.  These  programs  could  provide  prediction  of  performance 
for  a  given  condition  within  the  limitations  of  the  data. 
These  solutions  were  for  the  specific  boat  tested  and  a 
scaled  family  of  geometrically  similar  boats. 

In  the  early  1970 's  attempts  for  steady  state  solution 
from  mathematical  models  began  to  be  published.   One  pioneer 
in  the  field  has  been  John  S.  Letcher,  publishing  parts  of 
his  model  in  references  (1) ,  (2) ,  and  (3) .   Further  work  has 
been  done  by  Myers  (4)  and  Kerwin  (5) .   The  computational 
technique  for  solution  of  the  force  equilibrium  condition 
has  been  described  by  Hazen  (6) . 

By  parameterizing  forces  in  four  degrees  of  freedom,  the 
model  developed  in  this  paper  grows  from  these  previous  works. 
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1 . 3  Overview 

In  order  to  resolve  the  aerodynamic,  hydrodynamic  and 

hydrostatic  forces  acting  on  a  sailing  boat  a  coordinate 

system  and  origin  are  established.   Achieving  an  equilibrium 

of  forces  and  moments  requires  that  the  summation  of  all 

forces  and  moments  acting  either  through  the  origin  and  about 

the  origin  are  equal  to  zero.   To  facilitate  modeling, 

component  forces  which  act  independently  or  nearly  so  are 

identified  and  examined.   Once  the  forces  acting  on  the  system 

and  the  location  and  direction  through  which  they  act  are 

identified,  each  may  be  modeled  and  superimposed  to  give 

the  required  equilibrium  solution. 

1.3.1  Coordinate  System 

Figure  1  presents  a  coordinate  system  which  describes 

the  sailing  boats  velocity  and  attitude.   The  coordinate 

system  designated  by  the  capital  letters  X  ,  Y  ,  and  Z   are 

located  on  the  free  surface  with  X  and  Y   axis  in  the 

o      o 

horizontal  plane.   The  X  axis  lies  in  the  direction  of  the 
c  o 

boat's  motion  and  the  Y   axis  is  to  starboard  at  right  angles 

o 

to  this  velocity.   The  Z   axis  is  vertically  downward. 

2  o 

A  second  coordinate  system  is  fixed  in  the  boat  at  the 
forward  end  of  the  waterline  as  presented  in  Figure  2.   This 
coordinate  system  is  designated  by  the  small  letters  x,  y, 
and  z.   The  x  axis  is  oriented  along  the  fore  and  aft  center- 
line  on  the  boat's  waterline  with  the  y  axis  in  the  transverse 
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Figure   I  -  Coordinate  Systems 
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Figure    2 
Boat   Coordinate   System 
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direction  to  starboard.   The  z  axis  desends  orthogonally  to 
the  x-y  plane  in  the  direction  of  the  keel. 

The  angular  differences  between  the  two  coordinate 
systems  describe  the  attitude  of  the  boat  as  it  moves  through 
the  water.   The  positive  direction  for  these  angles  are  also 
presented  in  Figure  2.   Their  direction  was  selected  to  pro- 
vide positive  values  when  the  boat  is  in  equilibrium  under 
normal  sailing  conditions  with  a  weather  helm.   The  order  in 
which  these  angles  are  taken  is  important  to  the  final  atti- 
tude.  The  boat  is  first  rotated  about  the  Z   axis  through 

o  3 

the  angle  \p .      The  angle  <p  is  identical  to  the  boat's  leeway 

angle.   The  boat  is  then  rotated  about  the  X  axis,  an  angle  $ , 

which  is  the  boat's  heel  angle.   The  boat  is  then  rotated 

about  the  Y   axis,  an  angle  0  .   9  is  the  angle  by 

which  the  boat  pitches  down  at  the  bow  and  is  also  the  angle 

of  the  axis  about  which  the  boat  heels.   Finally  the  angle  a_ 

is  the  angle  the  rudder  is  deflected  from  the  boat's  centerline. 

A  sailing  boat  responding  to  the  forces  acting  upon  it 

and  conforming  hydrostatically  to  deformations  in  the  free 

surface  would  heave  in  the  positive  Z   direction  an  amount  s. 

*  o 

The  selection  of  these  coordinate  systems  allows  the 
modeling  of  forces  in  either  coordinate  system.   When  modeling 
velocity  related  forces,  such  as  wave  drag  and  lift,  the 
coordinate  system  fixed  in  the  water  is  more  attractive 
computationally.   On  the  other  hand,  hydrostatic  forces  are 
more  easily  computed  in  the  coordinate  system  fixed  in  the 
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boat.   In  the  final  solution  of  force  equilibrium,  forces 

and  moments  are  resolved  in  the  X  ,  Y  ,  and  Z   system  except 

o   o       o  ■*  * 

the  heeling  moment  which  is  resolved  about  the  x  axis. 

It  is  a  common  practice  in  naval  architecture  to  describe 
a  vessel's  dimensions  and  coefficients  in  a  positive  sense. 
With  few  exceptions  this  is  done  throughout  this  paper. 
The  positive  direction  of  forces  and  moments  is  defined  in 
Section  1.3.2. 

1.3.2  Equilibrium  of  Forces 

The  forces  and  moments  acting  upon  a  sailing  boat 

in  these  coordinate  systems  may  be  divided  into  six  components 

X  =  Forces  in  the  positive  X  direction 

r         o 

Y  =  Forces  in  the  positive  Y   direction 

^         o 

Z  =  Forces  in  the  positive  Z   direction 

e  o 

K  =  Moments  about  the  x  axis 

M  =  Moments  about  the  y  axis 

N  =  Moments  about  the  z   axis 

o 

The  positive  direction  of  the  moments  K,  M,  and  N  are  taken 
in  the  same  sense  as  the  angles  <f> ,    0 ,  ty      respectively  as 
indicated  in  Figure  2.   Under  normal  conditions  sail  aero- 
dynamic forces  will  create  positive  forces  and  moments  on 
the  boat. 

Through  out  the  remainder  of  this  paper  the  moments  K, 
M,  and  N  may  be  referred  to  interchangeably  as  forces  or 
moments . 
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The  resolution  of  forces  into  equilibrium  requires  that 
the  summation  of  all  the  forces  acting  in  the  six  directions 
be  equal  to  zero,  i.e.; 


EX.  =  EY.  =  EZ.  =  EK.  =  EM.  =  EN.  =  0 
111111 


where  "i"  designates  individual  aerodynamic  or  hydrodynamic 
forces  or  moments  acting  in  the  prescribed  direction.   An 
example  of  a  contributing  force  is  X  ,  the  sail  driving  force. 
This  is  one  of  a  number  of  forces  acting  along  the  X   axis. 

The  achievement  of  an  equilibrium  of  forces  in  the  six 
directions  requires  the  solution  of  six  non-linear  simultaneous 
equations.   To  obtain  a  unique  solution,  six  variables  are 
required.   A  selection  of  variables  for  the  six  degrees  of 
freedom  problem  is: 


(Ub,  if/,  s,  <(>,  9,   ar) 


The  only  unidentified  variable  is  U,  ,  which  is  the  boat's 
velocity  in  the  positive  X  direction.   Each  variable  contri- 
butes a  major  effect  to  a  corresponding  force.   The  boat's 
resistance  is  strongly  dependent  on  boat  speed,  U,  .   As  lee- 
way angle,  ip,  strongly  contributes  to  side  force,  so  does 
heel  angle,  <j>,  to  the  righting  moment  and  rudder  angle,  a  ,  to 
the  yaw  balance. 
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The  present  model  excludes  the  solution  in  heave  and 
pitch.   This  exclusion  of  two  degrees  of  freedom  was  not  due 
to  their  lack  of  importance,  but  the  lack  of  an  adequate 
parametric  model.   In  reality  the  axis  about  which  the  boat 
heels  may  be  pitched  as  severely  as  15°.   This  pitch  makes 
a  measurable  change  in  the  boat's  angle  of  incidence  when 
heeled  to  large  angles.   Further,  due  to  the  drop  in  dynamic 
pressure  as  the  water  passes  around  the  hull  a  certain  amount 
of  sinkage  will  occur  as  a  measure  of  the  hydrostatic 
reaction.   Besides  the  effect  on  form  drag,  this  dynamic 
pressure  and  resulting  sinkage  also  increases  the  boat's 
wetted  surface  and  viscous  resistance.   The  magnitude  of 
these  effects  is  small  enough  that  they  may  be  reasonably 
ignored.   In  the  generation  of  mathematical  models  for  the 
remaining  forces,  emperical  data  was  used  that  came  from 
experimental  tests  which  included  the  effects  of  both 
pitch  and  heave. 

The  present  model  fixes  the  pitch  angle,  9,  and  sinkage, 
s,  at  zero.   The  variables  are  reduced  to  the  following  four: 


(  Ub,  ip,  c|>,  ar) 


The  remaining  major  forces  acting  on  a  sailing  boat  are 
resolved  into  X,  Y,  K,  and  N  components  and  are  added  together 
to  give  the  following  force  equilibriums: 
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EX.  =  X   +  X.  +  X,  +  X   =0 

1     s     h    k     r 


SY.=Y   +Y.+Y   =0 
l    s    k    r 


ZK.=K   +  K.  +  K,  +  K   =0 
l     s     h    k     r 


ZN.=N    +  N    +N,+N   =0 
l    sx    sy    k    r 


The  individual  forces  are  described  as  follows 


X  =  Sail  driving  force 
s  3 


X,  =  Hull  resistance 
h 


X,  =  Hull  and  keel  induced  drag 


X   =  Rudder  induced  drag 


Y   =  Sail  side  force 
s 


Y,  =  Hull  and  keel  side  force 
k 


Y   =  Rudder  side  force 

r 


K   =  Sail  heeling  moment 


K,  =  Hull  righting  moment 


K,  =  Keel  heeling  moment 


K  =  Rudder  heeling  moment 
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N   =  Yawing  moment  from  sail  driving  force 

N   =  Yawing  moment  from  sail  side  force 
sy        J 

N,   =  Yawing  moment  from  keel  side  force 

N   =  Yawing  moment  from  rudder  side  force. 

Each  individual  force  is  defined  in  later  sections  as  a 
function  of  the  boat's  geometry,  a  true  wind  vector  and  the 
four  variables  U,  ,  i> ,    <j>,  and  a  .   The  solution  of  equations 
1  through  4  for  a  given  boat  geometry  and  wind  condition 
provides  values  for  boat  speed  and  attitude.   From  these 
values  the  hydrodynamic  performance  of  the  specific  sailing 
boat  may  then  be  evaluated. 
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II.   FORCE  AND  MOMENT  MODELS 

2 . 1  Hydrodynamic  Forces 

The  hydrodynamic  forces  acting  on  a  sailing  boat  may  be 
conveniently  divided  into  four  areas  for  modeling.   The  four 
areas  are  hull  forces,  keel  forces,  rudder  forces,  and  hydro- 
static forces.   Hull  forces  consist  of  resistance  from  wave 
making,  viscous  effects,  and  form  drag.   Keel  forces  are 
made  up  of  lift  on  the  keel  and  hull,  and  drag  induced  by 
this  lift.   Rudder  forces  in  the  same  way  consist  of  lift 
and  induced  drag  from  the  rudder.   Finally,  hydrostatic 
forces  are  caused  by  the  change  in  position  of  the  center 
of  buoyancy  in  the  hull.   Each  of  these  areas  is  presented 
individually  and  models  are  developed  for  their  contribution 
to  the  total  force  equilibrium. 
2.1.1  Hull  Resistance 

The  estimation  of  hull  resistance  is  a  problem  faced 
by  naval  architects  in  the  design  of  every  ocean  going 
vehicle.   The  difference,  in  predicting  a  sailing  boat's 
resistance,  is  in  the  determination  of  the  effects  changes  in 
attitude  have  on  these  calculations.   The  approach  taken  is 
to  determine  a  model  for  up-right  resistance  and  modify  this 
resistance  as  necessary  to  account  for  changes  in  the  boat's 
attitude . 
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Upright  Resistance 

The  technique  commonly  used  for  prediction  of  resistance 
from  tow  tank  data  is  based  on  Froude ' s  Hypothesis.   Basically, 
the  hypothesis  states  that  the  total  upright  resistance  may 
be  divided  into  two  components.   The  first  is  frictional  or 
viscous  resistance.   Frictional  resistance  is  dependent  on 
Reynold's  Number,  a  non-dimensional  parameter  defined  as: 


R  =  U.  L/y  „ 
n    b  '  rH20 


The  second  component  is  residual  resistance.   The  residual 
resistance  is  the  combination  of  drag  due  to  the  generation 
of  waves  and  the  effects  of  eddy  making  and  form  drag.   The 
effect  of  the  hull  generating  waves  is  a  gravitational  effect, 
largely  dependent  on  Froude  Number.  This  non-dimensional 
coefficient  is  defined  as: 


F   =  U, //gL 
n    b'  ^ 


Though  Froude  left  the  separation  of  the  elements  of 
resistance  in  these  two  components,  form  drag  and  eddy  making 
resistance  have  little  dependence  on  either  Froude  or  Reynolds 
Numbers  over  the  velocity  ranges  encountered  by  sailing 
boats.   For  this  reason  it  is  logical  to  further  separate 
residual  resistance  into  wave  resistance  and  form  resistance. 
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The  separate  components  of  upright  resistance  are  assumed 
independent  of  each  other  so  that  the  total  resistance  of  a 
boat  at  a  given  speed  may  be  estimated  by  superposition  of 
the  forces.   Total  upright  resistance  is  then: 

Upright  resistance  =  frictional  resistance 
+  wave  resistance  +  form  resistance. 

Each  component  of  upright  resistance  is  examined  in  detail 
and  mathematical  expressions  are  developed  to  quantify  their 
magnitude. 

Wave  Resistance 

The  mathematical  model  for  wave  resistance  is  based  on 
a  generalized  form  of  Michell's  Integral.   This  expression 
is  evaluated  for  a  base  hull  form  and  is  corrected  for 
parametric  deviation  from  this  parent  form.   This  expression 
is  emperically  adjusted  to  give  accurate  prediction  of  the 
towing  tank  results  of  Antiope  (3) . 

Michell  (7)  postulated  that  wave  resistance  for  thin 
ships  depends  upon  the  computation  of  pressure  changes,  6  , 
due  to  the  wave  pattern.   His  expression  for  this  resistance 
is 


=  !££_  /  (I2+J2) 


*w  "   „  2 


dA 


*Ub  X         A2-l 
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where 

+e   K(x)  ~  ,2       /TT2  , 

I  -     /    /  |Z     e"A   ^Z/Ub      cos  (iS*)    dx   dz  6 

-e    0  dX  IT 

D 

e      K(x)     a  ,2       /rT2  , 

J  =   /      /  *-*-  e        3        b        sin  (-^sr)    dx  dz  7 

-e    0  8x  U^ 

b 

Another  form  of  these  equations  was  derived  by  Weiblum  (8) 
A  general  expression  of  his  result  is 


2  2 

B  H     L 
RT  -  Pg — t —  F  (sr~ i    form  coefficients,  F  ) 
w       Li  j.  n 

c 


By  assuming  that  the  block  coefficient,  C,,  varies  very  little 
between  yacht  forms,  equation  8  may  be  written  as: 


3 

R.  =  —  [ y]   f  (m— /  form  coefficients,  F  ) 

W    Pg   (L/100)J       c  n 


The  number  of  experiments  conducted  to  validate  Michell's 
result  is  small.   Wigley  (9)  compared  model  results  to  theory 
and  observed  the  interference  humps  in  the  resistance  curve 
were  overpredicted  by  theory.   Further  comparison  between 
vessels  of  constant  draft  and  varying  beam  done  by  Wigley  (10), 
Weinblum  (11) ,  and  Mumford  (12)  indicate  for  slender  ships 
wave  resistance  is  approximately  proportional  to  the  square 
of  the  beam  as  estimated  by  Michell's  Integral,  but  as  the 
beam  to  length  increases  dependence  of  resistance  on  beam 
appears  to  be  less  than  quadratic.   It  is  expected  that  the 
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high  beam  to  length  ratios  encountered  in  sailing  boats  will 
strain  the  slender  ship  limits  of  Michell's  Integral. 

The  basic  form  of  the  wave  resistance  is  rewritten  as: 


R  =  —  [ o]nf  (L/T  ,  form  coefficients,  F  )       10 

w    p9   (L/100)3       c  n 


where  n  is  between  one  and  two.   The  actual  value  is 
determined  by  emperical  data  discussed  in  the  following 
sections.  «* 

The  identification  of  the  most  important  parameters 
affecting  wave  resistance  is  required.   Todd  (13)  did  a 
regression  analysis  of  over  two  hundred  destroyers  and  found 
the  following  coefficients  had  the  highest  correlation  between 
the  speed  length  ratios  of  0.8  and  1.25. 


(V/ (L/100)3,  C  ,  LCB,  L/Tc,  Fr) 


For  geometrically  similar  boats  we  know  directly  from 

Froude ' s  Hypothesis  that  the  wave  resistance  scales  is  the 

cube  of  the  characteristic  length.   The  general  form  of 

Michell's  Integral,  equation  10,  is  used  for  the  dependence 

on  displacement  to  length  ratio.   By  assuming  base  value  for 

prismatic  coefficient,  C   ;  longitudinal  position  of  the 

center  of  buoyancy,  LCB  ;  and  length  to  draft  ratio,  L/T   ; 

an  expression  for  the  wave  resistance  of  a  parent  form  may 

be  written: 
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T3         A 

R    =  f .  (Uu//L)  #-  [ S =-]n  11 


«P     X   b      pg   (L/100)3 


at  C   ,  LCB  ,  and  L/T 

po      o  CO 


To  predict  the  wave  resistance  of  a  sailing  boat,  the 
resistance  of  the  parent  boat  is  adjusted  by  correction 
factors  based  on  the  boat's  geometric  deviation  from  the 
parent.   The  wave  resistance  equation  is 


R  =  R    (1  +  C   correction  factor  +  LCB  correction 
w    wp       p 

factor  +  L/T   correction  factor)  12 


Each  component  of  this  equation  is  defined  and  discussed 
in  the  following  sections. 

Characteristic  Length  -  L 

The  characteristic  hydrodynamic  length  or  "effective 
length"  of  a  sailing  boat  is  difficult  to  quantify.   This 
has  been  a  difficulty  in  the  maintenance  of  all  offshore  rating 
rules  as  demonstrated  by  the  extremes  to  which  designers  have 
gone  to  make  rated  length  look  shorter  than  the  boat's  actual 
hydrodynamic  length.   For  the  purposes  of  this  model/  L  is 
determined  to  be  the  load  water  line,  LWL.   This  creates 
little  difficulty  when  examining  the  comparitive  performance 
of  similar  forms,  but  subjective  adjustment  would  be  required 
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when  comparing  plumb  bow,  transom  stern  boats  to  long  ended 
boats  with  shallow  buttocks  and  entrance  profiles.   To  some 
extent  the  latter  might  increase  in  length  with  Froude  Number 
as  the  boat  will  be  supported  more  by  her  overhangs. 
Parent  Froude  Number  Dependence  -  f ,  (U,//L) 

The  speed  length  ratio,  U,  //L",  was  adopted  instead  of 
the  dimensionally  correct  variable  F   =  U,//gL  due  to  the 
common  usage  and  the  comfortable  feel  most  yacht  designers 
have  with  speed  length  ratio.   The  function,  f, (U,//L)  adjusts 
the  wave  resistance  equation  to  accurately  predict  the 
resistance  of  the  yacht  Antiope .   The  function  was  determined 
by  inputting  Antiope ' s  geometry  into  the  model  and  selecting 
f ,  (U,  //L)    to  insure  proper  correlation  with  the  full  scale 
towing  test  results  reported  in  (3)  and  (14).   Values  of  the 
function  were  established  at  each  data  point  given  in  (14) 
and  a  spline  cubic  curve  was  fitted  to  give  the  curve  in 
Figure  3  for  n  =  2.0. 

Displacement  to  Length  Factor  -  [A/(L/100)  ] 
The  high  beam  to  length  ratios,  B/L,  encountered  in 
sailing  boats  requires  that  the  exponent,  n,  be  determined 
from  emperical  data.   Unfortunately,  the  data  base  from  tests 
of  actual  sailing  boats  is  small.   The  scale  effects  pointed 
out  in  reference  (15)  precludes  the  use  of  results  from  towing 
tests  of  small  scale  yacht  forms. 
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Figure    3 
Parent  Froude    Number  Dependence 
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Examination  of  Taylor  Standard  Series  (16)  at  prismatic 
coefficients  typical  in  yacht  design  indicate  that  a  value 
of  n  between  1.8  and  2.0  appears  reasonable.   Though  the 
prismatic  coefficients,  length  to  draft  and  displacement  to 
length  ratios  in  reference  (16)  are  typical  of  sailing  boats, 
the  block  coefficient,  C,  ,  is  high  and  not  characteristic  of 
yacht  forms . 

Figure  4  is  a  plot  of  two  curves  predicting  the  wave 
resistance  of  a  Standfast  43  sailboat  by  the  mathematical 
model  with  n  equal  one  and  two.   Also  plotted  are  the  data 
points  from  the  towing  tests  of  a  one  sixth  scale  model 
conducted  by  M.I.T.  and  Delft  University.   Figure  5  compares 
the  full  scale  test  (15)  of  the  one  ton  yacht  Bullet  to  two 
theoretical  curves.   The  form  drag  data  from  both  yachts  was 
taken  as  the  residuary  drag  at  a  speed  length  ratio  equal  to 
0.4  and  subtracted  from  the  residuary  resistance  versus  speed 
curve  to  obtain  the  wave  drag. 

Two  is  maintained  as  the  value  of  n  in  the  computerized 
model  appended  to  this  paper.   Though  comparisons  made  in 
Figures  4  and  5  indicate  that  n  =  2  may  slightly  overestimate 
the  effect  displacement  to  length  ratio  has  on  wave  making 
resistance,  the  limited  sample  size  available  for  correlation 
does  not  justify  adjustment  from  the  theoretical  value. 
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C   Correction  Factor 

The  prismatic  correction  factor  corrects  the  magnitude 

of  the  sailing  boat's  wave  resistance  for  deviations  away 

from  the  base  prismatic  coefficient,  C   .   The  base  prismatic 

po  c 

coefficient,  C   ,  was  selected  to  be  the  optimum  prismatic 

po  c  c 

coefficient,  C     ,  for  each  speed  length  ratio.   Taylor's 

Standard  Series  (16)  data  are  the  only  tests  available  of 

large  models  where  prismatic  coefficient  was  systematically 

varied.   Examination  of  these  data  reveals  a  definite  optimum 

prismatic  coefficient,  C    ,  ,  for  each  tabulated  speed  length 

ratio.   Figure  6  is  a  spline  cubic  fit  of  this  optimum 

prismatic  taken  at  the  displacement  to  length  ratio  equal  to 

250  and  plotted  versus  speed  length  ratio.   This  optimum 

prismatic  coefficient,  C    ,  ,  as  a  function  of  the  speed 
^  popt  c 

length  ratio  is  selected  as  the  base  prismatic  coefficient, 

C   ,  of  the  parent  equation.   This  is  done  so  that  deviations 

from  this  optimum  will  always  result  in  an  increase  in  wave 

making  resistance. 

From  Taylor's  Standard  Series,  reference  (16),  the 

percentage  increase  in  resistance  for  deviation  from  the 

optimum  prismatic  coefficient  was  tabulated  for  each  speed 

length  ratio.   From  this  emperical  data  the  expression  for 

the  C   correction  factor  is  generated. 
P 

C   correction  factor  =  f  2  (Ub//L)  1 100  (Cpopt'S*  ' 1  *  ?      13 
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Figure  7  is  a  plot  of  the  spline  cubic  fit  of  the  correction 
factor's  dependence  on  speed  length  ratio,  f „ (U,//L) .   A 
reasonable  definition  for  the  sailing  boat's  prismatic 
coefficient,  C  ,  must  be  established  to  insure  quantitative 
agreement  between  the  Taylor  Series,  reference  (16),  forms  and 
the  yacht  forms.   A  method  for  weighting  the  volume  and 
sectional  area  of  the  keel  in  the  calculation  of  the  longi- 
tudinal prismatic  coefficient  is  developed.   Equations  6  and 
7  of  Michell's  Integral  indicate  that  displacement  should  be 
weighted  exponentially  in  the  vertical  direction  by 

e-A  gz/ub. 

There  are  two  common  practices  in  calculating  the  C   of 
a  sailing  boat.   The  first  is  to  include  the  total  volume 
and  sectional  area  of  the  keel  into  the  calculation.   The 
inclusion  of  all  the  deeply  placed  sectional  area  of  the  keel 
would  overestimate  its  effect  on  the  wave  making  resistance 
when  compared  to  the  keel's  actual  effect  on  the  wave  train. 
On  the  other  hand,  the  second  method  is  to  treat  the  keel 
as  an  appendage  and  calculate  only  a  fairbody  or  canoe 
prismatic.   This  ignores  the  keel's  displacement  effect 
altogether.   Typical  values  of  C   for  the  same  yacht  calculated 

ir 

both  ways  are  0.52  and  0.58  with  and  without  the  keel 

included  respectively.   If  the  Taylor  Series,  reference  (16), 

optimum  is  to  be  used,  the  corresponding  prismatic  coefficient 

would  be  somewhere  in  between  these  two  values. 
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Figure   6  -  Optimum  Prismatic  Coefficient 
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Figure    7 
Prismatic    Correction    Factor 
Dependence    on   Froude    Number,      fp 
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One  would  expect  the  volume  of  the  keel  to  have  its 

largest  effect  when  the  sailing  boat  is  sitting  on  a  wave 

of  approximately  the  same  length  as  the  waterline.   At  this 

point  the  wave  trough  would  bring  the  keel  the  closest  to 

the  free  surface.   An  estimation  of  keel's  proper  effect  on 

prismatic  coefficient,  C  ,  may  be  made  by  weighting  the  keel's 

area  and  volume  by  e    nz   where  F   is  the  vertical  Froude 

2  nz 

number  evaluated  at  speed  length  ratio  equal  to  one.   The 
longitudinal  prismatic  coefficient,  C  ,  is  defined  as: 

V   +  V* 

C   =  — 14 

P    (A   +  A*)  L 
m    k 

where 


V*  =  V   e-32.2d/(1.689L)2  15 

k    k 


and 


A*  =  A   e-32.2d/(1.689L)2  16 

k    k 

V   and  V,  are  the  displaced  volume  of  the  canoe  body  and 
c      k  c 

keel  respectively.   A   and  A,  are  the  maximum  cross  sectional 
area  of  the  canoe  body  and  the  cross  sectional  area  of  the 
keel  at  that  station  respectively.   The  depth,  d,  is  measured 
from  the  static  waterline  to  the  center  of  the  keel's  displaced 
volume . 
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L/T  Correction  Factor 

Taylor  Standard  Series  (16)  tabulates  data  on  ships  at 
length  to  draft  ratios  equal  to  eight  and  sixteen.   Taylor 
Series  (16)  asserted  that  linear  interpolation  between  these 
ratios  was  adequate. 

The  base  length  to  draft  ratio,  L/T   ,  of  16  is 
selected  from  Taylor  Series  (16)  data  and  the  correction 
factor  is 


L/T   correction  factor  =  -f , (U.//L) (16-L/T  )  17 

C  j   D  C 


where  f3(U.//L)  is  a  spline  cubic  fit  of  the  speed  length 
effect  as  shown  by  Taylor  (16).   f  (U.//L)  is  plotted  in 
Figure  8. 

The  major  differences  between  the  ships  tested  for 
Taylor's  Standard  Series  (16)  and  yacht  forms  are  in  the 
shape  of  the  sections.   The  Taylor  form  typically  has  a 

maximum  section  coefficient/  C  ,  of  0.98  while  sailing  boats 

m  3 

range  between  0.60  and  0.70.   As  a  result,  the  vertical 

distribution  of  volume  between  the  two  forms  differ  to  a 

large  extent.   It  is  not  completely  accurate  to  assume 

that  vertical  shifts  in  the  volumetric  distribution  of  the 

two  forms  will  have  a  similar  effect  on  wave  making  resistance 

The  model  utilizes  the  best  data  available  to  date,  but  the 

need  exists  for  a  large  scale  yacht  series  examining  just 

such  a  systematic  parameter  variation. 
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L/Tc    Correction    Factor 
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LCB  Correction  Factor 

In  the  same  manner  as  the  C   correction  factor,  the  LCB 

P 

correction  factor  is  based  on  the  actual  LCB  and  an  optimum 
LCB  position.   A  summation  of  recommended  optimum  LCB  positions 
as  a  function  of  block  coefficient  and  prismatic  coefficient 
was  presented  by  Todd  (17) .   There  is  reasonable  agreement 
among  the  recommended  curves  with  less  than  1.5%  variation 
between  extreme  curves  in  the  area  of  interest.   Comparing 
the  LCB  curves  with  common  practice  in  yacht  design  a  curve 
attributed  to  Professor  Troost  in  (17)  was  adopted.   In 
modified  form  it  is 


LCB   .  =  62.5  -  17. 5C  18 

opt  p 


where  LCB     is  measured  in  percent  L  from  the  forward 

opt  e 

perpendicular . 

Examination  of  the  major  series  in  references  (17) ,  (19) , 
and  (20)  which  consider  LCB  as  a  variable,  none  adequately 
cover  the  region  in  which  sailing  boats  fall.   In  all  cases, 
these  series  cover  block  coefficients  higher  than  typical 
yachts.   Certain  trends  may  be  determined  though: 

1.  As  with  the  previous  factors,  there  is  a  peak 
effect  at  speed  length  ratio  equal  to  one. 

2.  As  block  coefficient  decreases,  the  penalty  for 

non-optimum  center  of  buoyancy  position  also 

decreases . 

The  above  series,  therefore,  place  an  upper  bound  on  the 
expected  penalty.  44 


Todd  (17)  indicates  for  Series  60  that  a  one  percent 
change  from  the  optimum  LCB  position  will  increase  wave 
making  resistance  by  approximately  one  percent.   Beyond  a  one 
percent  change  the  penalty  increases  sharply. 

The  assumed  factor  is 

LCB  correction  factor  =  sin  (90U,//L)  |  62  .  5-17  .5C  -LCB|-/100 

19 

Frictional  Resistance 

The  mathematical  model  for  frictional  resistance  is 
derived  from  procedures  adopted  by  the  Eighth  International 
Towing  Tank  Conference  (21),  ITTC  in  1957.   Traditionally 
frictional  or  viscous  resistance  is  defined  as: 


Rf  =  Cf.l/2p  Sw  V2b 


S      is  the  wetted  surface  of  the  sailing  boat  and  C,.  is  the 
w  3  f 

coefficient  of  frictional  resistance. 

The  ITTC  defined  the  coefficient  of  frictional 
resistance  as: 


Cf  =  0.075/(log10  Rn  -  2)2 


where  Reynold's  Number,  R  ,  is  the  only  parameter  affecting 
this  coefficient.   Again,  Reynold's  Number  is  defined  as: 
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R   =  U,  L/Y„  n 
n    b    Ho 


where  yu  „  is  the  kinematic  viscosity  of  water. 
H2° 

A  suitable  definition  of  the  Reynold's  Number  character- 
istic length,  L,  must  be  established.   The  ITTC  adopted  the 
use  of  the  load  waterline  length,  LWL,  as  the  characteristic 
length.   This  makes  sense  for  ships  with  rectangular  profiles, 
but  sailing  boats  characteristically  have  cut  away  profiles 
and  keel  and  rudder  chords  which  have  little  dependence  on 
the  length  of  the  waterline.   Davidson  (22)  realized  this 
problem  and  adopted  seventy  percent  of  the  waterline  length 
as  the  characteristic  length.   The  adoption  of  a  constant 
length  lacks  sensitivity  to  the  effects  differing  profiles 
have  on  the  characteristic  length.   To  make  the  mathematical 
model  sensitive  to  these  variations,  the  total  wetted  surface 
of  the  yacht  is  divided  into  three  areas,  each  with  its  own 
characteristic  length. 

The  three  areas  are  the  sailing  boat's  canoe  body,  the 

keel  and  the  rudder.   The  load  waterline  length,  LWL,  is  the 

characteristic  length  of  the  canoe  body  while  the  mean  chord 

of  the  keel,  c,  ,  and  rudder,  c  ,  are  the  characteristic 
k  r 

lengths  of  the  keel  and  rudder  respectively. 

The  coefficient  of  friction  for  each  area  is  calculated 
using  equation  21.   The  total  frictional  resistance  is  the 
sum  of  the  resistances  of  each  area  and  is  written  as: 
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R_    =    l/2p„    „    ul     (C„,S..,  +C.rtS1+C,r    S    ) 

f    '   H?0   b    ffb  fb   fk  k   fr  r 


The  above  coefficients  and  wetted  surfaces  are  defined  as: 
C.p.p,  =  fairbody  coefficient  of  friction 

Sf,   =  fairbody  wetted  surface 

Cr,   =  keel  coefficient  of  friction 
fk 

S,    =  keel  wetted  surface 
k 

C.e   =  rudder  coefficient  of  friction 
fr 

S    =  rudder  wetted  surface 

r 

The  ITTC  equation  for  the  coefficient  of  friction  is  an 
emperical  formulation  which  historically  has  either  slightly 
overpredicted  or  underpredicted  frictional  resistance  for 
certain  forms.   From  previous  experience  with  a  given  form, 
a  factor  may  be  added  to  give  satisfactory  correlation  with 
full  scale  results.   When  equation  22  is  applied  to  Antiope  (3), 
the  correlation  factor  is  unnecessary.   This  would  indicate 
that  the  yacht  Antiope  has  no  form  resistance  or  skin  roughness. 
This  is  unfounded  and  indicates  that  the  ITTC  formulation 
overpredicts  the  frictional  resistance  of  the  yacht  Antiope 
by  the  approximate  magnitude  of  its  own  form  resistance.   This 
effect  is  discussed  further  in  the  next  section. 
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Form  Resistance 

Form  resistance  is  a  combination  of  effects  which  show 
little  dependence  on  Froude  Number  and  Reynold's  Number  in 
the  speed  range  seen  by  sailing  boats.   Form  resistance,  R-  , 
is  defined  as :  • 


R*   =  C-   l/2pS„  U^  24 

fm    fm  '     H  fb  b 


The  form  resistance  coefficient,  C,  ,  is  a  constant  with 

fm 

speed  and  depends  only  on  the  boat's  geometry,  surface 
roughness,  and  appendages.   Form  resistance  is  made  up  of  eddy 
making  resistance,  roughness  effects,  and  pressure  drag  due 
to  separation. 

The  form  resistance  of  a  yacht  may  constitute  ten  percent 
of  the  boat's  total  resistance  at  low  speed  and  decreases  to 
less  than  two  percent  as  the  effects  of  wave  making  resistance 
increase.    Kirkman  and  Pedrick  (15)  discuss  the  difficulty  in 
modeling  form  resistance  from  analytic  and  emperical  sources. 
The  scale  effects  from  small  scale  model  tests  have  made  the 
measurement  of  this  component  difficult.   Further,  the 
selection  of  the  ITTC  friction  coefficient  formulation  includes 
an  arbitrary  amount  of  form  resistance  which  is  not  quantifi- 
able.  The  use  of  this  line  in  predicting  Antiope  (3)  resistance 
did  not  require  the  use  of  a  form  drag  coefficient  to  provide 
correlation. 
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The  present  mathematical  model  does  not  parameterize 
the  form  resistance  coefficient  but  allows  the  user  to 
input  a  value  to  improve  correlation  with  known  results. 
The  inclusion  of  this  coefficient  allows  future  parameteriza- 
tion as  techniques  for  its  prediction  become  available. 

Hull  Force  Formulation 

The  total  upright  resistance,  R  ,  is  the  superposition  of 
wave  resistance,  frictional  resistance,  and  form  resistance 
as  defined  by  equations  12,  23,  and  24. 


R  =  R  +  R.  +  R_  25 

u    w    f     fm 


The  change  in  hull  resistance  with  heel  and  leeway  angle  must 
be  considered.   For  small  angles  of  leeway,  ip ,  the  change  in 
resistance  is  considered  an  effect  of  lift  and  is  modeled  in 
the  section  on  keel  forces,  2.1.2. 

In  the  absence  of  lift,  the  change  in  upright  resistance 
due  to  heel,  <j> ,  is  caused  by  a  number  of  effects.   As  the 
sailing  boat  heels,  one  would  expect  not  only  changes  in  the 
boat's  symetrical  lines,  but  also  her  form.   Depending  on  the 
yacht  form,  the  effective  sailing  length  may  also  increase  or 
decrease.   Further,  as  the  boat  heels,  the  keel  moves  closer 
to  the  free  surface  and  in  effect  decreases  the  effective 
prismatic  coefficient  and  increases  resistance. 
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Examination  of  a  number  of  towing  tank  tests  conducted 

at  M.I.T.  show  that  the  increase  in  resistance  at  zero  side 

force  varies  depending  on  the  design.   This  variation  ranged 

from  an  increase  in  resistance  from  a  minimum  of  five  percent 

to  a  maximum  of  twenty-five  percent  at  twenty-five  degrees 

heel.   The  majority  of  the  well-designed  models  show  an 

increase  of  approximately  ten  percent  at  twenty-five  degrees 

heel.   Because  of  the  lack  of  test  date  available  which 

systematically  investigates  heel  effects,  the  increase  in 

resistance  with  heel  is  limited  to  an  estimate  of  what  good 

design  practice  has  demonstrated.   This  effect  is  emperically 

modeled  by  dividing  the  upright  resistance  by  the  cosine  of 

the  heel  angle. 

The  expression  for  the  hull  forces  in  the  X  direction 

o 

is  written  as: 


X.  =  R  /cos  d>  2  6 

h    u     Y 


where  R   is  the  upright  resistance  defined  by  equation  25. 

2.1.2  Keel  Forces 

The  mathematical  models  for  lift  forces  acting  on  the 
hull  and  keel  are  derived  from  slender  body  theory  as  applied 
to  fish-like  bodies  in  references  23,  24,  and  25.   Expressions 
are  defined  for  the  hull  and  keel  side  force,  Y, ;  induced  drag, 
Xj_;  yawing  moment,  N,  ;  and  the  heeling  moment,  K,  . 
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Hull  and  Keel  Side  Force 

Letcher  (3)  investigated  the  possibility  of  utilizing 
available  theories  for  predicting  a  sailing  boat's  side  force. 
He  compared  lifting  line  theory,  slender  body  theory  and  an 
adoption  of  slender-body  theory  for  fish-like  bodies  to  the 
towing  tank  test  results  of  the  yacht  Antiope .   The  theory 
presented  by  Newman  and  Wu  (24)  for  fish-like  bodies 
demonstrated  excellent  correlation  to  Letcher's  (3)  regression 
analysis  of  Antiope ' s  test  data.   An  adaptation  of  this  theory 
is  utilized  for  modeling  the  hydrodynamic  lifting  forces 
acting  on  the  keel  and  hull. 

The  side  force  acting  on  the  keel  and  hull  is  defined  as: 


Yk  =  -Cy*cos2* (1/2pH20  L2ub*'  27 

2 

where  the  cos  $    term  was  found  by  Letcher  (3)  to  satisfactorily 

quantify  the  effect  heel  has  on  the  side  force  coefficient, 
C  /  for  the  yacht  Antiope.   L  is  the  load  waterline  length. 

A  mathematical  model  for  side  force  coefficient  is 
derived  from  slender-body  theory  as  applied  to  fish-like  forms. 
Milgram  (23)  gives  an  expression  for  the  local  lift,  L(x),  on 
a  slender-body  as: 


L(x)  =  ip  |-  [b2(x)]V,U\     x   <  x  <  x,  28 

2MH^0  dx         b  b     o        b 
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where  b (x)  is  the  local  span,  x   is  the  horizontal  position  of 
the  maximum  span,  and  x,  is  the  horizontal  distance  from  the 

D 

maximum  span  to  the  forward  perpendicular.   Figure  9  shows 
Antiope's  profile  with  the  locations  of  x   and  x,.   V.  is 
the  component  of  velocity  y  direction. 

The  total  side  force  predicted  by  slender-body  theory  is 


y   =  r 

sb    J      L(x)  dx  29 

c 


For  the  numerical  computation  of  the  slender-body  side  force, 
the  yacht's  profile  is  tabulated  at  the  five  points  shown  in 
Figure  9  and  the  profile,  b (x) ,  is  approximated  by  straight 
lines  between  these  points.   Equation  29  is  evaluated 
numerically  to  give  the  slender-body  lift  on  the  hull  and 
keel. 

The  velocity  V,  may  be  approximated  for  small  leeway 
angle  by  the  linear  expression: 


Vb  =  V  30. 


The  coefficient  of  side  force,  C  ,  .  ,  based  on  slender-body 

y^sb  J 

theory  may   be  written   as: 


'    ,    u   =   Y    u/ip„    r,   L2    U?iL>  31 

yij;sb  sb/  2KH   0  by 


where   L   is   taken  as    the  waterline    length,    LWL. 
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Newman  and  Wu  (24)  developed  a  procedure  for  determining 
the  effect  an  axisymetric  body  of  radius  r   has  on  the  slender- 
body  lift.   From  reference  24,  Figure  10  is  a  plot  of 

2 
-L/irptLV,  b   versus  r  /b  .   At  r  /b   equal  to  zero,  the  curve 
/Kbbo         o/o       o   o   ^ 

gives  the  slender-body  lift  for  a  given  span.   As  r  /b 

increases  the  lift  is  attenuated  by  the  amount  shown  in  the 

curve.   This  curve  may  be  defined  as  a  fair  body  factor,  fft. 

The  curve  is  entered  with  the  ratio  r  /b   to  obtain  the  fair 

o/    o 

body  factor.   The  slender-body  coefficient  from  equation  31 
is  multiplied  by  the  fair  body  factor  to  compute  the  yacht's 
side  force  coefficient.   The  expression  for  the  side  force 
coefficient  is 


C    .    =   C    ,    .  (f„  )  32 

y\p  yijjsb      fb 


The  radius,  r  ,  is  measured  as  the  draft  of  the  canoe  body 

at  the  maximum  draft  measurement,  b  .   These  two  measurements 

o 

are  depicted  in  Figure  9. 

This  technique  for  the  solution  of  side  force  correlates 
well  with  the  full  scale  tests  of  the  yacht  Antiope  (3). 
Antiope  has  a  moderately  long  keel  compared  to  modern  sailboats 
with  their  fin  keels.   For  a  given  hydrodynamic  span,  as  the 
chord  of  a  long  keel  decreases,  the  lift  remains  essentially 
unchanged,  but  at  a  certain  point  a  further  decrease  in  the 
keel's  chord  can  no  longer  support  the  spanwise  elliptic 
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Figure    10  -  Lift   Attenuation 
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loading  assumed  by  slender-body  theory.   The  present  model  is 
not  sensitive  to  this  effect,  and  though  the  model  is  adequate 
for  modeling  changes  in  draft,  it  should  not  be  utilized  to 
optimize  keel  area. 

Keel  Induced  Drag 

Newman  in  reference  25  develops  the  equation  for  induced 

drag  on  fish-like  forms.   Figure  11  is  from  reference  25  and 

qives  a  curve  for  (^Y./D.)  as  a  function  of  the  ratio  r  /b  . 
y  k   l  o   o 

By  generating  a  spline  cubic  fit  to  this  curve,  induced  drag, 
X^.*  is  defined  as: 

xk  =  Yk^/(^Yk/Di)  33 

Keel  Yawing  Moment 

The  yawing  moment  due  to  the  side  force  on  the  keel  is 
defined  as : 


N,  =  x   .  Y.  34 

k    cpk  k 


The  center  of  effort  of  the  hydrodynamic  side  force,  x  .  ,  is 
determined  from  slender-body  theory.   Milgram  (2  3)  gives  the 
yawing  moment  about  x   in  Figure  9  as: 


x, 
Nsb  =  fD   xL(x)dx  35 

x 

o 
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Figure   11  -    Induced    Drag 
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where  L (x)  is  defined  by  equation  28.   Hence,  the  non-dimen- 
sional coefficient  for  yawing  moment  can  be  written  as: 


C  ,  u  =  N  ,  /^-pHo0  L3  uj>  36 

nipsb    sb  2   2      br 


and  the  distance  to  the  center  of  hydrodynamic  effort  from 
the  forward  perpendicular,  x  .  ,  as: 


x   .  =  (1.0  -  C  ,  ,/C  .  .  )L  37 

cpk  nipsb  yipsb 


Keel  Heeling  Moment 

The  heeling  moment  caused  by  the  side  force  on  the  keel 
is 


K,  =  -Y,  z   ,  co'sijj  38 


where  z   ,  is  the  distance  to  the  vertical  hydrodynamic  center 

of  effort  measured  from  the  waterline.   Since  the  leeway  angle, 

\p ,    is  small,  cosip  is  taken  as  unity. 

Slender-body  theory  tends  to  overpredict  the  location  of 

the  vertical  center  of  effort  when  compared  to  towing  tank 

data.   Letcher  (3)  reduced  Antiope  data  for  the  location  of 

the  center  of  effort  and  deduced  that  the  vertical  center  of 

effort  was  approximately  1.5  feet  below  the  waterline.   From 

this  a  reasonable  estimate  for  the  vertical  center  of  effort, 

z     is  made. 
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z   .  =  0.33b  39 

cpk        o 


where  b   is  the  maximum  draft  as  defined  in  Figure  9. 
o 

2.1.3  Rudder  Forces 

Mathematical  models  for  the  forces  acting  on  the 
rudder  are  derived  from  the  emperical  equation  for  spade 
rudders  presented  in  reference  (26)  .   Expressions  are  defined 
for  the  rudder  side  force,  Y  ;  rudder  drag,  X  ;  rudder  yawing 
moment,  N  ;  and  rudder  heeling  moment,  K  . 

Rudder  Side  Force 

Whicker  and  Fehlmer  (26)  conducted  a  series  of  wind  tunnel 
tests  on  low  and  medium  aspect  ratio  spade  rudders  and  pro- 
vided emperical  equations  which  accurately  model  the  forces 
and  moments  acting  on  a  large  family  of  spade  rudders.   The 
equations  are  sensitive  to  aspect  ratio,  a  ,  sweep  angle,  ft, 
and  taper  ratio,  t  . 

The  modeling  of  rudder  forces  in  relation  to  a  sailing 

boat  are  complicated  by  the  effect  of  the  downwash  from  the 

keel.   By  treating  the  keel  as  a  low  aspect  ratio  lifting 

surface,  the  downwash  can  be  assumed  to  deflect  the  flow 

along  the  centerline  of  the  boat.   Although  the  keel  will 

actually  deflect  the  flow  slightly  less  than  this,  the 

approximation  is  considered  reasonable.   The  model  therefore 

assumes  that  the  angle  of  attack  on  the  rudder  and  the  rudder's 

deflection  from  the  boat's  centerline,  a  ,  are  identical.   The 

effect  of  heel  on  rudder  performance  is  considered  negligable. 
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Figure  12  is  a  drawing  of  a  spade  rudder  including  the 
measurements  and  nomenclature  used  in  determining  the  rudder's 
performance.   The  geometrical  ratios  required  by  the  theory 
are  derived  from  these  measurements.   The  hydrodynamic  aspect 
ratio,  a  ,  is  defined  as: 


a   =  2b  /c  40 

r     r   r 


where  b   is  the  geometric  span  and  c   is  the  mean  chord  as 
defined  in  Figure  12.   Taper  ratio,  t  ,  is  defined  as: 


t  =  c  ./c  41 

r    rt  rr 


where  c   is  the  rudder  tip  chord,  and  c    is  the  rudder  root 

rt  r  rr 

chord  as  defined  in  Figure  12. 

Whicker  (26)  defines  the  coefficient  of  lift  as: 


dC  C,    a   .  « 


where  C,   is  the  cross  flow  drag  coefficient,  and  the  rudder 
angle,  a   is  measured  in  degrees. 

An  emperical  expression  from  reference  (27)  for  the  cross 
flow  drag  coefficient,  C,  ,  for  rudders  with  faired  tips  is 


C,   =  0.08  +  0.72t  43 

dc  r 
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direction  to  starboard.   The  z  axis  desends  orthogonally  to 
the  x-y  plane  in  the  direction  of  the  keel. 

The  angular  differences  between  the  two  coordinate 
systems  describe  the  attitude  of  the  boat  as  it  moves  through 
the  water.   The  positive  direction  for  these  angles  are  also 
presented  in  Figure  2.   Their  direction  was  selected  to  pro- 
vide positive  values  when  the  boat  is  in  equilibrium  under 
normal  sailing  conditions  with  a  weather  helm.   The  order  in 
which  these  angles  are  taken  is  important  to  the  final  atti- 
tude.  The  boat  is  first  rotated  about  the  Z   axis  through 
the  angle  \p .      The  angle  \p    is  identical  to  the  boat's  leeway 
angle.   The  boat  is  then  rotated  about  the  X  axis,  an  angle  <J> , 
which  is  the  boat's  heel  angle.   The  boat  is  then  rotated 
about  the  Y  axis,  an  angle  0  .   0  is  the  angle  by 
which  the  boat  pitches  down  at  the  bow  and  is  also  the  angle 
of  the  axis  about  which  the  boat  heels.   Finally  the  angle  a_. 
is  the  angle  the  rudder  is  deflected  from  the  boat's  centerline. 

A  sailing  boat  responding  to  the  forces  acting  upon  it 
and  conforming  hydrostatically  to  deformations  in  the  free 

surface  would  heave  in  the  positive  Z   direction  an  amount  s. 

*  o 

The  selection  of  these  coordinate  systems  allows  the 
modeling  of  forces  in  either  coordinate  system.   When  modeling 
velocity  related  forces,  such  as  wave  drag  and  lift,  the 
coordinate  system  fixed  in  the  water  is  more  attractive 
computationally.   On  the  other  hand,  hydrostatic  forces  are 
more  easily  computed  in  the  coordinate  system  fixed  in  the 

20 


rudder  stock 


hull  line 
s 


^~crr 


Figure    12  -    Rudder    Geometry 
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8CT/3a      is   defined    (26)    as: 

Li  r 


3CL         (0.9)  (2tt)    ar 


(57.3)  [cosp>  r  4      +4      +   1.8] 
cos   ft 


44 


where  ft  is  the  sweep  angle  as  defined  in  Figure  12. 

The  rudder's  drag  coefficient,  C-,  ,  is  defined  as: 


2  2 

C-,   =  C ,   +  C,   Aa  e  +  0.0166  C,  45 

dr    do    lr  ' "  r  lr 


where  e  is  an  efficiency  factor  established  at  0.9  (26). 
C,   is  the  rudder  section  form  drag  coefficient  and  if  not 
specified  C,   is  assumed  to  be 


CdQ  =  0.0065  46 


2 
The  0.0166  C,    term  is  an  emperical  increase  in  functional 

resistance  due  to  lift  recommended  by  Kerwin  (28) .   This  form 
drag  coefficient  is  from  the  NACA  0015  section  (26) .   The 
interference  drag  caused  by  the  interaction  between  the  fair- 
body  boundary  layer  and  the  rudder  root  chord  is  not  accounted 
for  in  the  present  model  and  should  be  examined  in  future 
studies . 

By  transforming  the  rudder's  lift  and  drag  coefficients 
from  the  boat's  coordinate  system  into  the  system  located  on 

the  horizontal  surface  of  the  water,  the  rudder  side  force 
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coefficient  is  written  as: 


C   =  (C,  cosip      -  C,  sinto)  coscf>  47 

yr     lr    r     dr    y  Y 


Hence,  the  rudder  side  force  is 


Y   =  -C   (~p„  _  A  U^)  4  8 

r     yr  2KK20   r  b 


where  A   is  the  rudder's  plan  form  area.   The  plan  form  area 
is  defined  as: 


A  =  b  c  49 

r    r  r 


Rudder  Drag 

The  rudder  drag  coefficient  is 


C   =  C,  sirni;   +  C-.  cosiJj  50 

xr     lr    y  dr    Y 


where  C.   and  C,   are  defined  by  equations  42  and  45 
respectively. 

The  equation  for  rudder  drag,  X  ,  is 


x  =  "C^^w  n  Ar-U^  51 

r     xr  2  H00  r  b 
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Rudder  Yawing  Moment 

To  determine  the  yawing  moment  caused  by  rudder  forces, 
the  location  of  the  center  of  effort  on  the  rudder,  with 
respect  to  the  rudder  stock,  is  calculated.   The  rudder 
stock  is  assumed  to  pass  through  the  quarter  point  on  the 
mean  chord  as  shown  in  Figure  12. 

The  moment  coefficient  about  the  quarter  chord  (26)  is: 


Cmo/4  =  f0-25  -  'SfVo'  <&'«  -0 

Li      Li  r   r 

2   a   l57.3; 


where 


8Cm        .    l.ll/a^+4+2 

( 2)      =  -  -  £ 53 

K%C      C  =0    2 
d^L  Ul  U    Z      4 (a  +2) 


The  normal  force  coefficient,  C  ,  is 

n 


C  =  C,  cosa   +  C,  cosa  54 

n    lr    r    dr    r 


The  location  of  the  center  of  effort  with  respect  to  the 
quarter  point  on  the  mean  chord  is 


C  —/A    c 

x  -/A  -  -  m^/4  56 

pc/4       C 
r  n 
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The  distance  from  the  forward  perpendicular  to  the  center  of 
effort  of  the  rudder,  xCDr'  is 


x    =  x  — , .    +  x  57 

cpr    pc/4    rs 


where  x   is  the  distance  from  the  forward  perpendicular  to 
the  rudder  stock. 

The  yawing  moment  due  to  rudder  forces  is 


N   =  x    Y  58 

r    cpr   r 


Rudder  Heeling  Moment 

From  reference  26,  the  distance  to  the  spanwise  center 

of  effort  from  the  root  chord,  z   ,  is 

ro 

(Clr(4/37r)    cosar  +  Cdr   sinar)br 

zro  =  C —  •  59 

n 

The  distance  from  the  waterline  to  the  center  of  effort  is 


z    =  z   +  z  60 

cpr    ro    rrc 


where  z    is  the  distance  from  the  waterline  to  the  root 
rrc 

chord  as  shown  in  Figure  12. 

The  rudder  heeling  moment,  K  ,  is 


Kr  "  zcPrClr  <KjO  ArUb>  61 
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2.1.4  Hydrostatic  Forces 

A  satisfactory  model  for  the  hydrostatic  righting 
moment  is  presented  by  Kerwin  in  reference  (5) .   In  modified 
form,  Kerwin 's  (5)  model  for  the  hydrostatic  righting  moment, 
N  ,  is  adopted. 


Nr  =  -0.00224  [A/(L/100)3]  L4(f>  (d-j+d^+c^F  ) 


+  2   sin*  62 

g  y 


where  L  is  the  waterline  length,  LWL,  and  d.  ,  d-,    and  d_  are 

non-dimensional  stability  coefficients  which  depend  on  hull 

form.   The  term,  z   sin*  accounts  for  the  distance  from  the 

g     y 

waterline  to  the  center  of  gravity. 

Each  coefficient  has  physical  significance  and  may  be 
input  by  the  model  user  or  generated  from  the  expressions 
given  below.   A  systematically  varied  yacht  series  based  on 
the  yacht,  Standfast  43,  is  presently  being  tested  jointly 
by  M.I.T.  and  Delft  University.   The  hydrostatic  data  from 
this  series  was  utilized  to  generate  parameterized  expressions 
for  each  of  the  coefficients. 

The  coefficient,  d..  ,  is  the  height  of  the  metacenter 
above  the  waterline  divided  by  the  load  waterline  length,  LWL. 
An  expression  for  d,  which  accurately  predicts  the  Standfast 
43  series  is: 

m    1119(BWL/L)3  .  0.335/(L/T  j  63 

1    [A/(L/100)J]  C 
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The  coefficient,  d_,  is  the  non-linear  effect  on  the  righting 
moment.   The  emperical  expression  for  d2  is 


d2  =  a    [0.0822 (BLW/L)  -  0.02162]  64 


The  coefficient,  d_,  quantifies  the  effect  the  wave 
train  has  on  the  stability  at  high  Froude  Numbers.   This 
effect  depends  on  the  boat's  beam  to  length  ratio,  BWL/L, 
and  the  fineness  of  her  bow  and  stern.   The  expression  for 
d^  is  an  emperical  fit  to  Standfast  43  data  and  is  only 
dependent  on  the  beam  to  length  ratio. 


d3  =  [7.8  -  56  ,350(BWL/L-0.28)2*75  (10  5)  BWL/L  >  0.28 
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and 


d.  =  7.8  x  10~5  BWL/L  -  0.28 

66 

The  location  of  the  center  of  gravity  below  the  water- 
line,  z  ,  may  be  specified.   If  not  specified,  z   is  assumed 
to  be 


z   =  0.025L.  67 

g 
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2 . 2  Aerodynamic  Forces  and  Moments 

The  aerodynamic  forces  are  divided  into  two  components , 
the  forces  acting  on  the  sailplan  and  its  associated  rigging, 
and  the  forces  acting  on  the  hull . 

2.2.1  Sail  Forces 

The  aerodynamic  forces  are  defined  in  terms  of  lift 
and  drag.   Lift  acts  at  right  angles  to  the  apparent  wind  while 
drag  acts  parallel  to  the  apparent  wind.   These  forces  are 
transformed  into  components  fixed  in  the  boat  and  water 
coordinate  system  defined  in  Figures  1  and  2.   Expressions 
are  derived  for  the  sail  driving  forces,  X  ;  the  sail  side 

S3 

force,  Y  ;  the  yawing  moment  from  the  sail  driving  force,  N   ; 

S  S  X 

the  yawing  moment  from  the  sail  side  force,  N   ;  and  the  sail 
heeling  moment,  K  . 

Wind  Triangle 

Computation  of  the  magnitude  and  direction  of  the 
apparent  wind  is  necessary  for  resolution  of  the  aerodynamic 
forces.   The  convention  used  in  this  model  was  adopted  by 
Hazen  (6)  and  defines: 


V   =  Upright  apparent  wind  velocity 


=  /(V  siny)2  +  (Ub+V  cosy)2  68 


and 
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3  =  Upright  apparent  wind  angle 


=  arctanf (Vtsiny/(Ub+VtcosY) ]  69 

where  V   is  the  true  wind  velocity  and  y    is  the  true  wind 
angle.  The  trigonometric  relationships  are  graphically  depicted 
in  Figure  13.   This  convention  differs  from  the  apparent  wind 
as  measured  by  instruments  on  a  sailing  boat  by  the  magnitude 
of  the  leeway  angle. 

In  the  case  of  the  sail  plan,  heel  angle  changes  the 
chord  wise  direction  and  magnitude  of  the  apparent  wind. 
Since  the  spanwise  component  of  the  wind  makes  a  negligible 
contribution  to  lift  and  drag,  the  apparent  wind  angle  and 
velocity  are  decreased  by  the  magnitude  of  the  spanwise 
velocity.   The  apparent  wind  velocity  is 


V  =  /V2  [  (coscbsinB  )  2  +  cos(3  ]  70 

a     au     y    u  u 


and  the  apparent  wind  angle  is 


3  =  arctan  [  (cos(})Sin3u)  /cos$u]  71 


The  resulting  lift  is  normal  to  the  mast. 
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Vj.cosT 


Figure    13  -  Apparent   Wind    Triangle 
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Lift  and  Drag 

A  drawing  of  a  representative  sail  plan  is  presented  in 
Figure  14  illustrating  the  measurements  utilized  by  the  sail 
force  model. 

Lift  is  defined  as: 


L.  -  c,  rip  .   V2  A  )  72 

s    Is  2  air   a   s 


where  A   is  the  total  geometric  area  of  the  sail  plan. 

The  lift  coefficient,  C,  ,  depends  on  such  parameters 

-L  o 

as  sail  aspect  ratio,  boom  height,  mast  diameter,  and 
apparent  wind  angle.   The  coefficient  of  lift  is  also  changed 
by  crew  optimization  of  trim  and  camber  depending  on  the  wind 
direction  and  magnitude.   The  parameterization  of  the 
coefficient  of  lift  was  not  done  for  this  model.   The  model 
user  may  select  the  lift  coefficient  or  the  model  will 
default  to  1.25  as  a  reasonable  value.   A  possible  source 
for  lift  coefficients  as  a  function  of  sail  plan  geometries 
is  given  by  Milgram  (29) . 

The  sail  drag  may  be  defined  as  the  summation  of  form 
drag,  frictional  drag,  and  induced  drag.   In  coefficient 
form,  the  expression  for  sail  drag  is 


C 

D   =  (C-   A   +  C.  A  )ip  .   V2  +  (-^4)  L2/^p  .   V2A     73 
s     fms  m    fs  s  2  air   a       2'   s' 2  air   a  s 

Lls 
71 


The  friction  drag  coefficient,  Cf    ,  has  the  assumed  (29)  value 
of  0.04.   The  form  drag  coefficient,  C~      ,  based  on  the  cross- 
sectional  area  of  the  mast,  A  ,  is  a  function  of  the  apparent 
wind  angle,  3,  but  to  date  little  data  is  available  on  which 
to  base  a  model.   For  this  reason  the  form  drag  coefficient 
may  be  specified  by  the  user  or  default  to  a  constant  equal 
to  0.4.   The  cross-sectional  area  of  the  mast  is  assumed  (29) 
to  be 


Am=  O.OKz   +  xe) 
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where  z   and  x   are  defined  in  Figure  14. 
p      e  3 

The  induced  drag  of  the  sail  plan  is  a  function  of  the 

sail  plan  geometry.   Milgram  (29)  recommends  the  use  of  the 

2 

ratio  C, .  /C,    as  a  measure  of  sail  efficiency.   Reference 
dis'  Is  J 

(29)  describes  a  study  of  systematically  varied  sail  plans 

using  lifting  surf  ace  theory.  From  this  study,  an  approximate 

2 
spline  cubic  fit  of  C , .  /C.    was  established  as  a  function 
e  dis   Is 

of  the  aspect  ratio  of  the  sail  plan,  a  ,  with  a  correction 
for  the  difference  between  the  aspect  ratio  of  the  main,  a  , 

and  the  aspect  ratio  of  the  fortriangle,  af.   The  ratio 

2 

C-,-  /C,    is  defined  as: 
dis   Is 


C,.  /C,  2  =  f(a  )  +  (a  -a.) (0.0625-0.00893  a  )  75 

dis'  Is       s      m   f  s 


for  a   >  ar 
m    f 
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or 


C , .  /C,    =  f(a  )  for  a   -  a,.  76 

dis'  Is       s  m    f 


The  aspect  ratio  of  the  sail,  a  ,  is  defined  as: 

s 


a  =  z?/A  78 

s    is 


where  z.  is  the  mast  height  as  shown  in  Figure  14.   The  aspect 
ratio  of  the  main  is  defined  as: 


a  =  z  /x  79 

m    P   e 


where  z   and  x  are  shown  in  Figure  14.   The  aspect  ratio  of 


the  fortriangle,  af,  is  defined  as 


a^  =  z./x.  80 

f    i'  j 


where  z.  and  x.  are  shown  in  Figure  14.   The  spline  cubic 

2 
function,  f (a  ),  is  given  in  Figure  15.   The  ratio,  C, . _/C,   , 

5  d  JL  S    J.  S 

is  generated  for  mast  head  rigs,  so  the  model  is  capable  of 

accepting  the  ratio  designated  by  the  user.   The  model  may  be 

used  for  other  geometries  by  selection  of  coefficients  from 

reference  (29) . 

The  above  lift  and  drag  coefficients  are  useful  until  the 

apparent  wind  hauls  aft  and  the  jib  and  shrouds  interfere  with 

the  main.   At  this  time  the  sail  plan  stalls  and  a  different 

73 


Figure   14  -  Sail   Plan   Measurements 
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Figure    15  -  Sail  Efficiency 
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model  must  be  used  between  the  apparent  wind  angle  of  90°  and 
180°.   Data  for  systematic  variations  to  sail  plans  in  the 
off-wind  condition  is  not  available.   The  following  simple 
form  for  lift  and  drag  when  the  apparent  wind  angle,  8,  is 
90°  or  greater  is  adopted. 


L   =  L  cos  (6  -  90°)  for  6  -  90°  81 

so    s 


and  the  off-wind  drag  is 


D    =  D   +  [(1.2/2)p  .   A  V2  -  D  )sin(6-90°) 
so    s        '   H air  s  a    s     M 


for  3  -  90°  82 

If  the  off-wind  lift  and  drag  coefficients  are  known  for  a 
given  point  of  sail,  they  may  be  specified  by  the  user.   When 
data  becomes  available  on  reaching  the  off-wind  coefficients 
as  a  function  of  the  apparent  wind,  the  function  may  be 
incorporated  into  the  model. 

2.2.2   Hull  Forces 

The  wind  forces  on  the  hull  during  normal  sailing 
conditions  may  account  for  only  one  or  two  percent  in  direct 
resistance  and  ten  percent  of  the  side  force,  but  as  the  wind 
increases  and  sail  is  shortened,  the  hull  forces  will  eventually 
dominate  the  system.   Very  few  experiments  have  been  conducted 
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investigating  the  effect  of  wind  on  ship  f orms .   Hughs  (30) 
conducted  tests  in  water  of  the  topsides  of  typical  tanker, 
liner,  and  cargo  ship  profiles.   To  date  no  tests  have  been 
conducted  on  the  sailboat  form.   Hughs  (30)  did  test  the 
above  water  shape  of  bare  cargo  hulls  and  his  emperical 
results  are  used  to  approximate  these  forces  on  a  sailing 
boat.   The  magnitude  of  the  aerodynamic  hull  force,  F,  ,  is 


Fu  =  0.001428  V2  (Ansin2e  +0.3A^cos26  )/cos(A-8)        83 
h  a   1     u     t     u 


where  A,  and  A  are  the  longitudinal  and  transverse  projected 
area  of  the  hull  respectively.   The  angle,  X,  is  the  angle 
between  the  centerline  of  the  boat  and  direction  from  which 
the  hull  force  acts.   The  angle  X  is  taken  from  the  experimental 
results  on  hulls  reported  in  reference  (30) .   A  spline  cubic 
was  fit  through  this  data  and  is  presented  in  Figure  16. 

One  effect  which  cannot  be  ignored  when  examining  aero- 
dynamic forces  on  the  hull  is  the  vertical  wind  velocity 
profile.   The  hull  of  a  sailing  boat  sees  considerably  less 
wind  velocity  than  the  sail  plan.   Curry  (31)  measured  the 
wind  velocity  at  different  heights  and  plotted  a  profile. 
By  assuming  a  height  of  two  feet  above  the  water  surface  as 
representative  of  an  adequate  mean  for  sailing  boats,  the 
true  wind  velocity  on  the  hull  was  estimated  from  Curry's 
wind  profile  to  be  0.55  times  the  true  wind,  V  .   The  wind 


77 


180 


0    20  40  60  60  100  120 140  160  180 


Pu 


Figure  16  -  Hull  Force  Angle 
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triangle  is  calculated  with  this  lower  true  wind  velocity 
and  the  hull  force  is  calculated  with  the  resulting  apparent 
wind  velocity  and  angle. 

Calculation  of  the  yawing  moment  requires  the  location 
of  the  center  of  pressure  of  the  aerodynamic  hull  force.   The 
distance  from  the  forward  perpendicular  to  the  center  of 
pressure  divided  by  the  load  waterline  length,  was  formulated 
from  data  in  reference  (30) . 


x   ,/L  =  0.2  +  0.00311  6  84 

cplr  u 


The  yawing  moment,  N,  ,  is 


Nh=  FhL(xoph/L)sinX  85 


Resolution  of  Aerodynamic  Forces  and  Moments 
The  total  aerodynamic  driving  force,  X  ,  is 


X  =  L  sin3  -  D  cos3  -  F.  cosX  86 

s    s        s        n 


and  the  side  force  is 


Y  =  (L  cos3  +  D  sin3)cos<}>  +  F,  sinX  87 
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The  heeling  moment  is 


K  =  Z      (L  cosR  +  D  sinft)  +  2.0  F.  sinA  88 

s    cpys   s        s  h 


where  z     is  the  distance  from  the  vertical  center  of  effort 
cpys 

of  the  sail  side  force  to  the  waterline.   An  expression 

estimating  the  location  of  the  vertical  center  of  effort, 

z     #  from  data  in  reference  (29)  is 
cpys' 


z     =  0.43  z.  +  z„  89 

cpys         i    fb 


where  z.  is  the  mast  height  and  zf,  is  the  freeboard  measured 
from  the  waterline  to  the  base  of  the  mast  as  shown  in 
Figure  14. 

The  magnitude  of  the  leeway  angle  is  ignored  in 
resolving  yawing  moments.   The  yawing  moment  due  to  aero- 
dynamic side  forces,  N   ,  is 


N   =  x     (L  cos3  +  D  sin$)  cosd)  +  N,  90 

sy    cps    s    M    s    M     y  h 


where  x    is  the  distance  from  the  horizontal  center  of 
cps 

effort  of  the  sail  plan  to  the  forward  perpendicular  and  is 
defined  (29)  as: 


x    =  x  +  0.02  (x.  +  x  )  91 

cps    m         j  e 
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where  x   is  the  distance  from  the  forward  perpendicular  to  the 

mast,  x.  is  the  distance  from  the  jib  tack  to  the  mast,  and 

3 

x  is  the  length  of  the  main  boom  as  shown  in  Figure  14. 
The  yawing  moment  due  to  the  driving  force  is 


N   =  z      X  sine})  92 

sx    cpxs   s 


where  z     is  the  distance  from  the  vertical  center  of  the 
cpxs 

driving  force  to  the  center  of  resistance.   z     is  estimated 

CPXb 

by 


z     =  0.48  z.  +  z.,  93 

cpxs         l    fb 
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III.   GEOMETRY  SYNTHESIS 


3.1  Parameter  Generation 


The  computerized  version  of  the  mathematical  model  has 
the  capability  to  generate  dimensions  of  a  satisfactory  design 
from  whicht parametric  trade-offs  may  be  made.   The  user  of  the 
program  may  specify  as  many  or  as  few  of  the  dimensions  as  he 
desires.   This  section  presents  the  equations  which  determine 
the  unspecified  dimensions. 

The  equations  do  not  generate  an  optimum  design.   The 
majority  of  the  expressions  are  estimates  by  the  author  of 
the  dimensions  of  a  typical  ocean  racer/cruiser  and  no  other 
justification  is  made  for  their  formulation. 

If  the  waterline  length  is  unspecified,  the  length  is 
set  at: 

LWL  =30.0  feet  94 


3 

The  default  displacement  to  length  ratio,  A/(L/100)  ,  is 


A/(L/100)3  =  225  95 


and  the  prismatic  coefficient,  C  ,  is 


C   =  0.53  96 

P 
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The  beam  to  length  ratio,  BWL/LWL,  is  defined  by  a 
linear  fit  through  the  mean  of  a  large  sample  of  boats 
presently  measured  under  the  International  Offshore  Rule. 
The  beam  to  length  ratio  is  defined  by: 

BWL/LWL  =  0.3633  -  0.002018  LWL  97 

By  assuming  that  the  base  boat  has  a  constant  block 
coefficient/  C,  ,  the  length  to  fairbody  draft  ratio  may  be 
written  as : 


LWL/T   =  10975  (BWL/LWL) / [A/ (L/100) 3]  98 

c 


As  a  percentage  of  the  waterline  length,  the  longitudinal 
position  of  the  center  of  buoyancy,  LCB,  is 

LCB  =53.0  99 

The  fluid  properties  are  established  at  the  following  values 

kinematic  viscosity  of  water: 

Yu  n  =  1.2791  x  10~5  100 

H2° 


density  of  water: 

p    =  1.99  101 
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density  of  air: 


p  .   =  0.00238  102 

air 


The  mean  chord  of  the  keel,  c,  ,  is  estimated  by: 


c,  =0.23  LWL  103 

k 


The  form  resistance  and  correlation  coefficient,  C,.  ,  is 

rm 

assumed  to  be: 


C   =  0.0  104 

fm 


The  estimate  of  fairbody  wetted  surface  gives  a 
satisfactory  prediction  of  the  Standfast  43  Series  and 
predicts  Antiope ' s  fairbody  wetted  surface  to  within  1.5%. 
For  most  modern  yacht  forms,  the  expression  should  give  a 
satisfactory  estimate  for  preliminary  design  purposes.   The 
expression  for  fairbody  wetted  surface,  S-.  ,  is 

S-.  -  (LWL)2{1.7/(LWL/T  )  +3  .  5xl0~5  [A/  (L/100)  3]  (LWL/T  )  } 

L  L)  C  C 


{0.0712 [ (BWL/T  )-4.0]  +  1.0}  105 


The  sail  area,  A  ,  is  established  at 

s 


A   =  (LWL)2  106 

s 
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The  dimensions  of  the  sailplan,  as  described  in  Figure 
14,  are  defined  by  assuming  the  aspect  ratio  of  the  fortri- 
angle,  a~,    and  the  aspect  ratio  of  the  main,  a  ,  are  equal 
to  three.   The  foot  of  the  jib  is  assumed  to  be  160%  of  the 
base  of  the  fortriangle.   It  follows  that  the  mast  height  is 


z.  =  /A  /0.4044  107 

l     s 


and  the  base  of  the  fortriangle  is 


z.  =  Z./3.0  108 

i  l 


The  main  hoist  length,  z  ,  is  assumed  to  be 


z   =  z./l.l  109 

P    i 


and  the  main  boom  then  is 


z   =  z  /3.0  110 

e    p' 


The  horizontal  positions  of  the  profile  measurements  as 
defined  in  Figure  9  are  as  follows : 


x  =0.0  111 

o 


x1  =  0.1  LWL  112 
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x2  =  0.2  LWL  113 


x3  =  0.3  LWL  114 


x.  =  0.55  LWL  115 

b 


The  maximum  draft,  b  ,  is  assumed  to  be  the  base  draft  as 
defined  by  the  International  Offshore  Rule  (33)  and  is 
calculated  by: 


b   =  0.146  LWL  +2.0  116 

o 


At  the  other  previously  defined  horizontal  positions,  the 
draft  measurements  are 


b,  =  0.98  b  117 

1         o 


b2  =  0.95  LWL/(LWL/T  )  118 


b-  =  0.7  LWL/(LWL/T  )  119 


b,  =  0.0  120 

b 


The  ratio,  r  /b  ,  as  defined  in  Figure  9,  is 
'  o'    o 


r  /b   =  b   LWL/(LWL/T  )  121 

o'    o    o    '     '  c 
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The  sailplan  is  positioned  in  the  boat  so  that  the 

center  of  effort  of  the  sailplan,  x    ,  coincides  with  the 

c  cps 

center  of  effort  of  the  keel,  x  ,  ,  as  defined  in  equation  37. 

The  horizontal  distance  from  the  forward  perpendicular  to  the 

sail's  center  of  effort,  x    ,  is  defined  as: 

cps 


x    =  x   .  122 

cps    cpk 


The  forward  overhang,  FOH,  is  not  determined  by  asthetics, 
but  coincides  with  the  placement  of  the  jib  tack.   This  point 
is  calculated  by  the  expression: 


FOH  =0.02  (x.  +x)+x.  -x  123 

j     e     j  cps 


The  freeboard,  z,.,  ,  is  derived  from  the  base  freeboard 

fb 

as  defined  by  the  International  Offshore  Rule  (33)  with  a 
constant  added  to  account  for  the  height  of  the  cabin  top 
as  shown  in  Figure  14.   The  freeboard  is  defined  as: 


z„     =  0.057  LWL  +2.0  124 

fb 


An  estimation  of  the  projected  lateral  area  of  the  hull, 


A1#  is 


A,  -  (2.1  LWL  +  FOH)zfb/2.0  125 
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and  the  project  transverse  area  of  the  hull,  A  ,  is  estimated 
by  the  equation: 


A   =  1.1  LWL(BWL/LWL) zfb  126 


The  wetted  surface  of  the  keel,  S,  ,  is  estimated  by: 


S,  =  2.1c,   [b   -  LWL/(LWL/T  )]  127 

K.  JC     O  C 


The  dimensions  of  the  rudder  are  determined  by  the 
following  equations.   The  distance  from  the  forward  perpen- 
dicular to  the  rudder  stock,  x   ,  is 


x    =  0.95  LWL  128 

rs 


and  the  distance  from  the  waterline  to  the  root  chord  is 


z   =0.2  LWL/(LWL/T  )  129 


The  root  chord  of  the  rudder  is  defined  by: 


c    =  0.075  LWL  130 

rr 


and  by  assuming  a  taper  ratio,  the  tip  chord  is  given  by: 


c  .  =  0.667  c  131 

rt  rr 
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The  mean  chord   follows   as: 


c      =    (c        +   c   J/2.0  132 

r  rr  rt 


An  estimate  of  the  span  of  the  rudder  is  given  by  the 
expression, 


b  =  0.85(b   -  z  )  133 

r         or 


and  the  sweep  angle,  U,    is  set  to  zero, 


Q   =  0.0  134 


The  rudder  area,  A  ,  is 


A  =  b  c  135 

r    r  r 


and  the  rudder's  wetted  surface  is  estimated  by 


S   =  2.1  A  136 

r        r 


89 


IV.   PROCEDURE 

4. 1  Equilibrium  Computation 

To  find  the  equilibrium  solution  of  the  forces  and 
moments  for  a  given  wind  direction  and  speed,  equations  1 
through  4  must  be  solved.   The  equations  are  rewritten  as: 


Rn  =  X   +  X.  +  X.  +  X   =0  137 

1    s    h    k    r 


R2  -  Ys  +  Yk  +  Yr  =  °  138 


R,  =  K   +  K.  +  K.  +  K   =0  139 

3    s    h    k    r 


R,  =  N    +N    +  N,  +  N   =  0  140 

4     sx    sy    k    r 


The  solution  of  these  equations  require  a  technique  for 
solving  non-linear  equations.   The  technique  utilized  is 
Newton-Raphson  Iteration  described  in  references  (34)  and  (6) 

R  is  an  "error"  vector  defined  by  the  components  R. 
through  R., 


R  =  [R1#  R2,  R3,  R4]  141 


where  the  vector,  R,  is  equal  to  zero  when  the  sailing  boat  is 

in  force  and  moment  equilibrium.   Due  to  the  iterative  nature 

of  the  Newton-Raphson  technique,  the  value  of  the  vector,  R, 

is  minimized  to  an  acceptable  tolerance. 
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P  is  defined  as  the  vector  of  assumed  values  of  the 
four  independent  variables, 


P  =  [Ub,  \\j,    <f>,  ar] 


142 


and  the  first  partial  derivative  of  the  component  "errors" 
are  defined  as: 


11 


*21 


31 


'41 


8R1 
8P^' 


8R2 

8P1' 
9R4 


12 


22 


'32 


42 


3RX 

3P^' 


^2 
3P2' 


8R3 
3P2' 


13 


'23 


'33 


43 


3R1 
9PJ' 

^2 
3P3' 

8R3 

W3' 

9R4 
3P3' 


3R. 
14  =  8P~ 


24 


3R. 
3P" 


3R. 
34    3P~ 


'44 


3P 
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These  partial  derivatives  are  computed  by  taking  first 
differences. 

By  defining  <$"  as  a  vector  of  required  incremental 
changes  to  the  independent  variables,  a  matrix  of  linearized 
equations  for  each  "error",  R,  through  R4  is  established. 


sll  S12  S13  S14 

S21  S22  S23  S24 

s31  S32  S33  S34 

S41  S42  S43  S44 


t             \ 

-Rl 

-R2 

"R3 

"R4 

V.              J 

144 
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Therefore: 


[6r  52,  53,  54] 


Sll  S12  s13  s14 
S21  s22  S23  S24 
S31   S32   S33   S34 


S41   s42   S43   S44 


'-1 

-Rl 

"R2 

"R3 

j 

,            4J 
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The  new  values  of  the  independent  variables  becomes: 


P1(i+1)  =  P1(i)  +  6± 


146 


P2(i+1)  =  P2(i)  +  62 


147 


P3(i+D  =  P3(i)  +  <53 


148 


P4(i+1)  =  P4(i)  +  64 


149 


Convergence  is  tested  by  examining  the  magnitude  of  the 
component  "errors",  and  the  vector,  P,  is  updated  until  the 
desired  tolerance  is  satisfied. 
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V.    RESULTS 

Figures  17  and  18  are  performance  polar  plots  of  two 
yachts.   Figure  17  shows  the  theoretical  boat  velocity  of  the 
yacht  Antiope  as  predicted  by  the  mathematical  model.   The 
velocity  is  plotted  as  a  function  of  the  true  wind  angle,  y, 
and  the  three  wind  velocities  5,  7.5,  and  10  knots.   For  this 
prediction  the  sail  area  was  assumed  to  be  400  square  feet 
and  the  center  of  gravity  is  0.61  feet  below  the  waterline. 
The  sailplan  geometry  and  placement  was  generated  by  the  model 
and  the  keel  fixed  rudder  was  modeled  by  placing  the  program's 
spade  rudder  at  the  trailing  edge  of  the  keel  and  setting  the 
wetted  surface  of  the  rudder  equal  to  zero.   Since  the 
correlation  of  many  of  the  model's  forces  are  derived  from 
data  of  the  full  scale  testing  of  Antiope  (3),  her  performance 
polar  is  included  in  the  results  even  though  the  model  is  not 
designed  for  accurate  prediction  of  rudder  forces  on  keel 
mounted  rudders . 

Figure  18  gives  the  performance  polar  of  the  model 
generated  sailing  boat  for  the  same  three  wind  velocities. 
The  geometry  of  each  boat  is  fully  described  in  the  computer 
printout  listed  in  Section  8.3.   The  two  example  output 
listings  in  Section  8.3  are  the  source  for  the  plots  in 
Figures  17  and  18. 
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Figure    17-  Anfiope  Performance  Polar 
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Default  Boat  Peformance   Polar 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

Two  important  criteria  should  be  considered  when 
evaluating  the  capabilities  of  a  system  model.   The  first  is 
the  verification  of  the  accuracy  of  the  predicted  performance, 
and  the  second  is  an  assessment  of  the  capabilities  and 
limitations  a  model  has  for  predicting  changes  in  performance 
due  to  variation  in  the  controlling  geometric  parameters. 

To  determine  the  accuracy  of  the  model  in  predicting 
total  performance,  the  performance  data  of  a  known  yacht  is 
required.   Unfortunately,  the  measurement  of  the  performance 
of  an  actual  yacht  under  sail  has  never  been  accomplished. 
To  provide  a  basis  for  correlation,  the  four  independent 
variables,  boat  speed,  leeway  angle,  heel  angle,  and  rudder 
angle  would  have  to  be  measured  on  a  full  scale  yacht  along 
with  true  or  apparent  wind  velocity  and  direction.   This  data 
is  not  available  to  date. 

Other  programs  exist  for  resolving  tow  tank  data  and 
assumed  sail  coefficients  into  performance  polars.   These 
programs  might  be  used  in  the  future  as  a  basis  for  comparison 
if  large  or  full  scale  towing  tank  data  is  used. 

The  use  of  other  performance  prediction  programs  for 
correlation  is  not  useful  unless  another  standard  exists  with 
which  to  measure  their  capabilities.   If  two  performance 
programs  agree,  the  agreement  may  be  that  both  will  give  the 
same  wrong  answer. 
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The  polar  plots  in  Figures  17  and  18  are  considered 
reasonable  only  because  they  correlate  subjectively  with  the 
author's  own  sailing  experience.   Though  each  component  force 
is  correlated  against  either  emperical  or  theoretical  data, 
a  means  of  examining  the  capability  of  the  mathematical  model 
for  the  total  system  should  be  developed  in  the  future. 

The  capability  of  a  mathematical  model  to  respond 
correctly  to  changes  in  boat  geometry  and  therefore  enable 
parametric  tradeoffs  requires  physical  accuracy  of  the 
individual  force  models.   When  the  individual  force  models 
were  presented,  the  limitation  and  capabilities  of  each  was 
discussed.   The  limitations  and  capabilities  of  each  model 
are  summerized  below. 

Hull  Resistance 

Hull  resistance  is  calculated  by  superposition  of 
frictional  resistance,  wave  resistance  and  form  resistance. 

The  frictional  resistance  coefficient  is  based  on  the 
ITTC  line  and  is  applied  to  the  fairbody,  the  keel  and  the 
rudder.   This  line  includes  an  unspecified  amount  of  form 
resistance.   Future  model  refinement  should  include  a 
frictional  coefficient  which  does  not  include  form  resistance 
and  a  mathematical  model  that  parameterizes  form  resistance. 
The  form  resistance  coefficient  is  not  parameterized  in  the 
present  model. 
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Wave  resistance  is  estimated  by  a  combination  of  the 
results  of  Michell's  theory,  Taylor  Standard  Series  tests 
and  full  scale  tests  of  the  yacht  Antiope.   The  resulting 
semi-emperical  mathematical  model  for  wave  resistance  needs 
to  be  compared  to  the  results  of  a  series  of  full  scale  yacht 
wave  resistance  tests. 

The  mathematical  model  for  'hull  resistance  is  considered 
adequate  for  predicting  calm  water  resistance.   To  predict 
a  boat's  performance  in  rough  water,  a  mathematical  model 
predicting  the  added  resistance  due  to  waves  should  be 
developed. 

The  change  in  resistance  due  to  heel  angle  is  also  not 
parameterized  and  was  modeled  as  an  average  of  what  good 
designs  can  expect.   The  investigation  and  incorporation  of 
heel  effects  would  be  worthwhile. 

The  present  resistance  model  is  considered  dependable 
for  parametric  variation  in  displacement  to  length  ratio 
and  wetted  surface. 

Keel  Forces 

The  hydrodynamic  side  force  on  the  hull  and  keel  are 
derived  from  slender-body  theory  modified  for  fish-like  forms 
(24) .   This  model  provides  accurate  prediction  of  forms  fixed 
with  moderate  keels  such  as  found  on  the  yacht  Antiope.   A 
method  for  determining  decrease  in  lift  as  the  keel  chord 
decreases  at  a  given  span  should  be  developed.   The  present 
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method  adequately  models  changes  in  performance  due  to  draft 
but  is  insensitive  to  reduction  in  keel  area  due  to  chord 
length  changes. 

The  interference  effects  between  the  hull  boundary 
layer  and  the  root  chord  of  the  keel  should  be  the  subject 
of  further  investigation. 

Rudder  Forces 

The  rudder  forces  are  modeled  from  the  free  stream  tests 
of  low  and  medium  aspect  ratio  lifting  surfaces  (26) .   The 
model  is  adequate  for  yachts  with  spade  rudders,  but  many 
yachts  are  being  designed  with  skeg  rudders.   A  model  should 
be  developed  for  rudder  forces  as  a  function  of  the  percentage 
flap  on  the  skeg  rudder.   The  fairbody  boundary  layer  inter- 
ference effects  on  the  rudder  have  been  ignored  in  the 
present  model. 

The  mathematical  model  is  expected  to  satisfactorily 
model  changes  in  the  spade  rudder  geometry. 

Hydrostatic  Forces 

Hydrostatic  forces  are  modeled  using  the  concept  of  the 
metacentric  height  and  are  adjusted  for  Froude  Number  effects. 
The  model  is  considered  sensitive  to  geometric  changes  though 
further  verification  of  the  Froude  Number  effects  is  desirable. 

To  accurately  model  the  righting  moment  on  smaller  yachts, 

a  function  representing  the  effect  of  crew  weight  on  righting 

moment  should  be  added. 
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Sail  Forces 

The  sail  force  mathematical  model  is  based  on  lifting 
surface  theory  presented  in  reference  (29)  .   The  model  is 
adequate  for  close  hauled  prediction  but  is  not  parameterized 
for  reaching  and  running  due  to  lack  of  data.   Work  is  needed 
to  determine  force  coefficients  as  a  function  of  the  apparent 
wind  angle  and  sail  plan  geometry  for  off  wind  conditions. 
Further  work  needs  to  be  accomplished  to  parameterize  the 
form  drag  coefficient  as  a  function  of  rig  geometry  and 
apparent  wind  angle. 

The  utilization  of  the  model  should  be  limited  to  light 
and  moderate  wind  conditions  until  the  optimization  variables 
for  reefing  and  flattening  are  included. 

The  effect  of  the  vertical  wind  gradient  on  heel  has 
not  been  included  and  should  be  added  in  the  future.   When 
added  resistance  in  waves  is  investigated,  it  will  also  be 
important  to  study  the  effect  large  pitching  motions  have  on 
sail  forces  and  their  relationship  to  the  added  resistance 

in  waves. 

i 

The  limitations  of  the  sail  force  model  in  its  present 
form  preclude  the  accurate  prediction  in  off  wind  and  heavy 
air  conditions. 

Aerodynamic  Hull  Forces 

The  wind  forces  acting  on  the  hull  are  modeled  from 
data  on  cargo  ship  hulls  (30) .   Though  the  model  is  derived 
from  a  different  hull  form,  the  prediction  of  these  forces 
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is  considered  reasonable.   Tests  on  yacht  hulls  would  help 
improve  the  correlation  for  this  mathematical  model,  but 
the  effect  on  the  final  prediction  should  be  minor. 

In  Conclusion 

The  present  form  of  the  mathematical  models  is  useful 
for  preliminary  design  in  determining  the  quantitative  effect 
changes  in  the  major  geometrical  parameters  have  on  a  boat's 
performance.   These  parameters  include  displacement  to 
length  ratio,  wetted  surface,  beam  to  length  ratio,  total 
draft,  rudder  geometry,  position  of  the  center  of  gravity, 
sail  plan  geometry  for  on  wind  conditions,  and  sail  area. 

The  development  of  the  model  has  been  useful  in  cata- 
loguing available  theoretical  and  emperical  sources  for 
modeling  forces  acting  on  a  sailing  boat  and  has  made  clear 
the  necessary  improvements  required  in  the  science  of  sailing 
boats  to  improve  prediction  capability.   For  the  naval 
architect,  this  first  generation  system  model  for  sailing 
boat  performance  may  be  used  as  a  file  to  which  improvements 
may  be  made  as  a  greater  understanding  of  the  science 
develops. 
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VIII.  APPENDICES 

8 . 1  Computer  Input 

The  required  data  cards  are  divided  into  two  data 
sets  described  below: 

Card  Set  Number  One 

Card  set  number  one  consists  of  the  input  points 
necessary  to  generate  the  spine  cubic  curves  required  by  the 
force  models.   There  are  twenty-four  cards  punched  in  FORMAT 
(8F10.0).   This  data  set  is  an  integral  part  of  the  program 
and  does  not  change  from  run  to  run.   Table  1  is  a  listing 
of  the  required  data  cards. 

Card  Set  Number  Two 

Card  set  number  two  is  repeated  for  each  boat 

evaluated  by  the  mathematical  model.   The  cards  control 

the  program  options  and  describe  the  geometry  of  a  given 

sailing  boat.   A  listing  of  the  card  set  used  to  describe 

the  yacht  Antiope  is  given  in  Table  2.   Table  3  is  an 

example  of  the  required  data  set  where  all  parameters 

are  defaulted. 

Card  Set  #2 

card  #1  IFAM,  NGAM,  NVT;  FORMAT  (8110) ;  no  default. 

col  10  If  IFAM  =  1  forces  and  moments  are  printed. 

If  IFAM  =  0  no  forces  and  moments  are 
printed. 

col  11  -  20         NGAM  =  number  of  true  wind  angles  desired. 

col  21  -  30         NVT  =  number  of  true  wind  velocities  desired 

105 


card  #2  ITMAX,  IPRINT,  N;  FORMAT  (10X,  13,  17X,  II, 

19X,  13) ;  no  default. 

col  11  -  13         ITMAX  =  maximum  number  of  iterations 

allowed  per  solution. 

col  31  If  IPRINT  =  1,  the  independent  variables, 

apparent  wind  data,  and  forces  and  moments 
will  be  printed  out  each  iteration. 

If  IPRINT  =  0,  only  the  convergent  or 
final  iteration  will  be  printed. 

col  53  N  =  4,  the  number  of  degrees  of  freedom 

card  #3  EPS1,  EPS2;  FORMAT  (10X,  E7.1,  13X, 

E7.1);  no  default. 

col  11  -  17         EPS1  =  required  tolerance  in  force 

equilibrium. 

col  31  -  37         EPS2  =  step  size  taken  for  computing 

first  differences. 

card  #4  (XOLD(I),  1=1,4);  FORMAT  (20X,  5F10.3); 

no  default. 

col  21  -  30         XOLD(l)  =  initial  assumed  value  of  boat 

velocity,  (feet  per  second) 

col  31  -  40         XOLD(2)  =  initial  assumed  heel  angle, 

(degrees) 

col  41  -  50         XOLD(3)  =  initial  assumed  leeway  angle, 

(degrees) 

col  51  -  60         XOLD(4)  =  initial  assumed  rudder  angle, 

(degrees) 

card  #5  (KVT(I),I  =  1,NVT);  FORMAT  (8F10.5); 

or  more  no  default. 

col  1-80  KVT(I)  =  desired  true  wind  velocities 

(knots) 

card  #6  OGAMMA(I),  I  =  1,NGMA);  FORMAT  (8F10.5), 

or  more  no  default 

col  1-80  DGAMMA(I)  =  desired  true  wind  angles 

(degrees) 
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card  #7  C(2),  0(7),  0(8),  C(9),  C(10);  FORMAT 

(8F10.5);  default  =  -1.0 

col  1-10  C(2)  =  waterline  length  (feet 

col  11  -  20  C(7)  =  displacement  to  length  ratio 

col  21  -  30  C(8)  =  prismatic  coefficient 

col  31  -  40  C(9)  =  length  to  draft  ratio 

col  41  -  50  C(10)  =  distance  from  FP  to  LCB  (%LWL) 

card  #8  C(ll),  FORMAT  (E14.7),  default  =  -1.0 

col  1-14  C(ll)  =  kinematic  viscosity  of  water 

card  #9  (0(1),  I  =  12,18);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(12)  =  density  of  water 

col  11  -  20  0(13)  =  mean  keel  chord  (feet) 

col  21  -  30  0(14)  =  mean  rudder  chord  (feet) 

col  31  -  40  0(15)  =  form  and  correlation  coefficient 

col  41  -  50  0(16)  =  fairbody  wetted  surface 

col  51  -  60  0(17)  =  keel  wetted  surface 

col  61  -  70  0(18)  =  rudder  wetted  surface 

card  #10  0(1),  0(3),  0(4),  0(5),  0(6);  FORMAT 

(8F10.5);  default  =  -1.0 

col  1-10  0(1)  =  beam  to  length  ratio 

col  11  -  20  0(3)  =  sail  area 

col  21  -  30  0(4)  =  heeling  moment  arm 

col  31  -  40  0(5)  =  pitching  moment  arm 

col  41  -  50  0(6)  =  distance  from  FP  to  rudder  stock 
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card  #11  (C(I),I  =  19,26);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(19)  =  close  hauled  lift  coefficient 

col  11  -  20  C(20)  =  rig  form  drag  coefficient 

col  21  -  30  C(21)  =  cross  sectional  area  of  rig 

col  31  -  40  C(22)  =  sail  frictional  coefficient 

2 

col  41  -  50  C(23)  =  induced  drag  coef/lift  coef 

col  51  -  60  C(24)  =  lateral  area  of  the  hull 

col  61  -  70  C(25)  =  transverse  area  of  the  hull 

col  71  -  80  C(26)  =  distance  from  CE  sail  to  FP 

card  #12  (C(I),I  =  27,30);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(27)  =  off  wind  drag  coefficient 

col  11  -  20  C(28)  =  off  wind  lift  coefficient 

col  21  -  30  C(29)  =  distance  for  CE  hull  to  FP 

col  31  -  40  C(30)  =  density  of  air 

card  #13  (C(I),  I  =  31,38);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(31)  =  xQ 

col  11  -  20  C(32)  =  x± 

col  21  -  30  C(33)  =  x2 

col  31  -  40  C(34)  =  x3 

col  41  -  50  C(35)  =  xfa 

col  51  -  60  C(36)  =  bQ 

col  61  -  70  C(37)  =  b1 


col  71  -  80         C(38)  =  b2 
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card  #14  0(39),  C(40);  FORMAT  (8F10.5);  default=-1.0 

col  1-10  C(39)  =  b3 

col  11  -  20         C(40)  ■  b, 

b 

card  #15  (0(1),  I  =  41,43);  FORMAT  (8F10.0); 

default  =  -1.0 

col  1-10  0(41)  =  draft  of  the  canoe  body/total  draft 

col  11  -  20         0(42)  =  hydrodynamic  lift  coefficient 

col  21  -  30         0(43)  =  \pX   lift/drag 

card  #16  (C(I),  I  =  45,52);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  0(45)  =  distance  from  LWL  to  rudder 

root  chord 

col  11  -  20  0(46)  =  rudder  root  chord 

col  21  -  30  0(47)  =  rudder  tip  chord 

col  31  -  40  0(48)  =  rudder  span 

col  41  -  50  0(49)  =  rudder  sweep  angle 

col  51  -  60  C(50)  =  rudder  aspect  ratio 

col  61  -  70  0(51)  =  rudder  mean  chord 

col  71  -  80         0(52)  =  rudder  3C,/3a 

V       r 

card  #17  (0(1),  I  =  53,56),  0(61);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  0(53)  =  rudder  cross  flow  drag  coefficient 

col  11  -  20  0(54)  =  rudder  3(^/3^ 

col  21  -  30  0(55)  =  rudder  form  drag  coefficient 

col  31  -  40  0(56)  =  rudder  area 

col  41  -  50  0(61)  =  distance  from  CE  keel  to  LWL 
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card  #18  (C(I),  I  =  57,60);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(57)  =  metacentric  height/LWL 

col  11  -  20         C(58)  =  Froude  Number  righting  moment 

coefficient 

col  21  -  30         C(59)  =  non-linear  hydrostatic  coefficient 

col  31  -  40         C(60)  =  distance  from  LWL  to  CG 

card  #19  (C(I),  I  =  62,69);  FORMAT  (8F10.5); 

default  =  -1.0 

col  1-10  C(62)  =  mast  height 

col  11  -  20  C(63)  =  base  of  the  fortriangle 

col  21  -  30  C(64)  =  main  hoist 

col  31  -  40  C(65)  =  main  boom  length 

col  41  -  50  C(66)  =  aspect  ratio  of  the  fortriangle 

col  51  -  60  C(67)  =  aspect  ratio  of  the  main 

col  61  -  71  C(68)  =  aspect  ratio  of  the  sail  area 

col  71  -  80  C(69)  =  distance  from  mast  aft  to  sail  CE 

card  #20  C(70),  C(71);  FORMAT  (8F10.5);  default=-1.0 

col  1-10  C(70)  =  forward  overhang 

col  11  -  20  C(71)  =  freeboard 

card  #21  KEND,  FORMAT  (II) ,  no  default 

col  1  If  KEND  =  1,  indicates  this  is  final 

data  set 

If  KEND  =  0,  indicates  other  data  sets 
to  follow 
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8.2  Computer  Listing 

The  main  program  and  all  necessary  subroutines 
utilized  to  compute  the  mathematical  model  are  listed 
in  this  section. 
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8.3  Example  Output 

Two  example  computer  outputs  are  listed  in  this  section 
The  first  is  the  yacht  Antiope  and  the  second  is  the  model 
generated  sailing  boat  which  is  produced  if  all  geometric 
parameters  are  defaulted. 
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ITH1X   * 

10 

IPRIMT  - 

0 

H       - 

4 

BPS1    = 

1.00-10 

EPS2   ■  , 

1.0E-02 

XOLD(1)  .  ..X0LD(    4) 

1.000000E   01  0.000000E-01  0.000000E-01  0.000000E-01 

3OL0TI0SS    ABE   COMPUTED    FOB    THE    FOLLOWING    TBOE    IIHD    SPEEDS    (KHOTS) 

5.00000        7.50000      10.00000 

SOLUTIOHS    ABE   COHPOTED    FOB    THE   FOLLOHING    TBOE    UIHD    ANGLES    FOB    EACH   IIHD   SPEED 

180.00000    160.00000    140.00000    120.00000    100.00000      80.00000      60.00000      45.00000       35.00000 
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LBL=  24.30000 

DISPLACESENT    TO    LEHGTH=  180.39990 

PRISBATIC-  0.53000 

L/H=  12.15000 

LCB-  53.50000 

KINEHATIC    VISCOSITI=         0. 1279100E-04 

DENSITY    OP    BATEB-  1.99000 

BEAN    KEEL   CH0BD=  7.12500 

BEAR    RUDDER    CHORD=  0.00000 

BO0GB8ESS    AND    COHHELATIOH    FACTOB=  0.00000 

FAIBBODT    BETTED    SUB?AC£=         118.61000 

KEEL    BETTED    S0RFACE=  41.38200 

BODDER    WETTED    SURFACE*  0.00000 

TOTAL    BETTED    SURFACE*  159.99200 

BBL  /LBL=  0.23900 

SAIL    ABEA  =  400.00000 

HEELING    BOBENT   AHfl=  16.27361 

PITCHING    BOBENT    ABB*  17.84613 

DISTANCE    PROB    FP'TO    BOCDEB    STOCK=  16.00000 

CLOSE    HAOLED    LIPT    COEF*  1.25000 

BIG   POBB    DRAG   COEF*  0.40000 

CROSS   SECTIONAL    ABEA   OP    BIG=  19.60001 

SAIL    PRICTIONAL   C0EP«  0.04000 

INDUCED    DRAG   COEP   /    LIPT  C0EP»*2=  0.13030 

LATERAL    AREA   OP    HOLL=  100.00000 

TRANSVERSE   AREA    OP    HULL*  27.00000 

DISTANCE    BETBEEN    CE    SAIL    AND    PP-  11.50000 

OPP    BIND    DRAG   COEP*  0.00000 
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OFF    BIND    LIFT  COEP=  0.00000 

DISTANCE    BETWEEN    CE  HULL   AMD    PP=  11.35764 

DENSITY    Of   AIB  =  0.00238 

HORIZONTAL    DISTANCE  PBOH    BAXIBOH    DBAFT   OP    PROFILE    BEASUBEBENTS 

0.00000         5.75000        8.51000  14.67000       16.91000 

PBOPILE    BEASOREBBBTS 
4.65000      '4.62000         2.12000         0.78000         0.00000 

DHAPT    OP    CANOE   BODT   /  TOTAL    DBAFT*  0.35000 

HYDBODYNAHIC    LIPT    C0EP=  0.09805 

LABDA   X    LIFT   /   DRAG-  2.06498 

DISTANCE    FBOB    LWL    TO   VEBTICAL   HYDBO    CE=  1.54845 

•**    BDDDEB   SECTION    ••* 

DISTANCE    FBOB    L«L    TO    BOOT  CHOBD*  0.40000 

BOOT   CHOBD=  1.82250 

TIP   CHOBD=  1.21561 

RUDDER    SPAN*  3.61250 

SHEEP    AHGL2=  0.00000 

ASPECT    BAIIO*  4.75625 

BEAN   CHOBD=  1.51905 

DCL    /  DALPHA*  0.06745 

CBOSS   FLOH    DBAG    COEF-  0.56024 

DCB   /   DCL=  0.21407 

FOBS    DRAG   COEF  =  0.00650 

RUDDER    A8EA  =  5.48758 

•**    RIGHTING    BOBENT    SECTION    •** 

BETACENTRIC    HEIGHT   /  LHL=  0.05711 

FROODE    NOBBER    SIGHTING    BOBENT   COEP=         0. 7800000E-04 
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MOB    LINEAB   HYDROSTATIC    COEP=      -0.  1 127451E-03 

DISTANCE    BELOB    LWL   TO    CG  =  0.60750 

***   SAIL    SECTION    *** 

I»"         31.15027  J»  10.U83U2 

P*  28.5911*  E»  9.53038 

AB   FOR    T3I=  3.00000  AB   HiIN=  3.00000 

AB   SAIL   PLAB=  2.U7280 

CE  SAIL    PLAN    APT    BA3T=  0.40028 

POBHARD   0VERHAHG=  -0.61630  ^ 

fREEBOARD*  2.75000 
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K?T»  5.00000      KBOTS 

KVA=  2.48617       KBOTS 

SUCCESSFUL    CONVERGENCE 

ITER         ■         S 

BOAT    SPEED=  2.S1383      KNOTS 

LEEWAY    ASGLE=  -0.00390       DEGREES 


DGAfl(1A  =         180.00000      DEGREES 
DBETTA=  179.99980      DEGREES 


HEEL   ABGLE=  0.00128      DEGREES 

RODDER    ABGLE=  0. 0  1191  DEGREES 


KTT=  5. 00000      KNOTS 

KV»=  2.60820       KBOTS 

SOCCESSPOL   CONVERGENCE 
ITEH        *         2 

BOAT   SPEED=  2.72703      KBOTS 

LEEHAr    ANGLES  0.02088      DEGREES 


DGAHHA=         160.00000      DEGREES 
DBETTA-  139.03030      DEGREES 


HEEL   ANGLE=  0.02578      DEGREES 

RODDER    ASGLE=  0. 0  1 149 DEGREES 


BO   COBTEBGEBCE 

KVT»  5.00000       KBOTS 

KVA=  3.24689       KBOTS 

SOCCESSPOL  COBTERGEBCE 

ITBB         »      11 

BOAT   SPEEDS  3.36880      KBOTS 

LEEWAY    ANGLES  0.14458      DEGREES 


DGAHHA"         140.00000      DEGREES 
DBETTA=  98. 17024      DEGREES 


HEEL    ANGLES  0.40965      DEGRESS 

RODDER    ANGLE=  0. 00311DEGREBS 


K7T=  5.00000       KNOTS 

KVA=  4.67690       KBOTS 

SOCCESSPOL   COBVEBGEBCE 


DGAHflA=         120.00000      DEGREES 
DBETTA*  67.79747      DEGREES 
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ITEH        •        4 

BOAT   SPEED= 
LEEWAY    ANGLE= 


4.26724      KNOTS 

0.41392       DEGBEES 


HEEL   ABGLE=  2.74890       DEGBEES 

RUDDER    ANGLE*  0. 19046DEGBEES 


KTT=  5.00000      KBOTS 

KYA=  6.28702      KIOTS 

SUCCESSFUL   CONVEBGENCE 
ITBB         ■         4 

BOAT   SPEED=  4.77718      KNOTS 

LEEWAY    ABGLE=  0.66322       DEGBEES 


DGAflflA=         100.00000      DEGBEES 
DBETTA*  51.55513      DEGBEES 


HEEL    ABGLE=  6.77472      DEGBEES 

BODDEB    AMGLE=  0. 66225DEGBEES 


KVTs  5.00000       KNOTS 

KVA=  7.53630       KBOTS 

SUCCESSFUL   CONVEBGENCE 
ITEB         »         3 

BOAT   SPEED=  4.83670      KNOTS 

LEEWAY    ASGLE=  0.97314       DEGBEES 


DGABHA=  80.00000       DEGBEES 

DBETTA"  40.79662      DEGBEES 


HEEL   ABGLE=  10.99982      DEGBEES 

BUDDEB    ANGLE'  1.  19417DEGBESS 


K?T»  5.00000      KNOTS 

KVA=  8.15836       KNOTS 

SUCCESSFUL  CONVEBGENCE 
ITEB         -         3 

BOAT   SPEEDS  4.41289      KNOTS 

LEEWAY    ANGLES  1.54980      DEGREES 


DGABBA=  60.00000      DEGBEES 

DBETTA*  32.05679      DEGBEES 


HEEL   ANGLES  13.76510      DEGREES 

BUDDEB    ANGLES  1.  49892DEGREES 
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KVT=»  5.00000      MOTS 

KV1=»  7.97554      KHOTS 

SUCCESSFUL    COHVEBGEHCE 

ITEB         ■        4 

BOAT    SPEED=  3.61350      KHOTS 

LEEHAT    ASGLE=  2.09007       DEGBEES 


DGABI1A=  45.00000      DEGREES 

DBETTA=  26.31445      DEGBEES 


HEEL    AHGLE=  13.57446       DEGBEES 

BODDEB    AHGLE=  1. 52272DEGBEES 


KVT*  5.00000       KHOTS 

KVA=  7.44743       KHOTS 

SOCCESSFOL   COHVEBGEHCE 
ITEB         -         4 

BOAT   SPEED=  2.77712      KHOTS 

LEEWAI    AHGLE=  3.95441       DEGBEES 


DGAflBA-  35.00000      DEGBEES 

DEETTA=  22.64893      DEGBEES 


HEEL    AHGLE=  12.01359      DEGBEES 

BODDEB    AHGLE=  1. 51425DEGBBES 


KVT=  7.50000       KBOTS 

K7*=  3.76861       KHOTS 

SOCCESSFOL   COHVEBGEHCE 
ITEB         *        4 

BOAT   SPEED=  3.73182      KHOTS 

LEEHAT    AHGLE=  -0.00395       DEGBEES 


DGABHA*         180.00000      DEGBEES 
DEETTA*         179.99980      DEGBEES 


HEEL    ABGLE=  0.00292       DEGBEES 

BODDEB    A8GLE=  0. 01209DEGBEES 


K7T=  7.50000       KHOTS 

KVA=  3.97397       KHOTS 

SOCCESSFOL   COHVEBGEHCE 


DGAHBA=         160.00000      DEGBEES 
DEETTA-         139.79790      DEGBEES 
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ITEB         »         2 

BOAT   SPE2D= 
LEE8AY    AHGLE= 


4.01511      KNOTS 

0.01904       DEGREES 


HEEL    AHGLE=  0.05447       DEGBEES 

BDODEB    ANGLE=  0.  01699DEGBEES 


KVT=  7.50000       KNOTS 

R?A=  4.91665      KMOTS 

SUCCESSFUL   CONVEBGEHCE 
ITEB         »         4 

BOAT   SPBED=  4.77972      KMOTS 

LEEWAY    AHGLE=  0.13080      DEGREES 


DGAHBA=         140.00000      DEGBEES 
DBETTA*         101.32590      DEGBEES 


HEEL    AHGLE=  0.86860      DEGBEES 

BODDEB    AMGLE=  0. 05872DEGBEES 


KVT=  7.50000       KMOTS 

KVA=  6.75196       KNOTS 

SUCCESSFUL    COMVEBGENCE 
ITEB         *         4 

BOAT   SPEED=  5.59371      KMOTS 

LEEWAY    AHGLE=  0.28889      DEGBEES 


DGAHHA=  120.00000      DEGBEES 

DBETTA=  74.14812      DEGBEES 


HEEL    AHGLE=  4.78607      DEGBEES 

BODDEB    ANGLE=  0. 56792DEGBEES 


KVf=  7.50000      KMOTS 

KVA=  8.77144       KNOTS 

SOCCESSFOL  CONVERGENCE 
ITER         »         4 

BOAT   SPEED=  6.03347      KMOTS 

LEEWAY    ANGLE=  0.35957       DEGBEES 


DGA.1BI=         100.00000      DEGBEES 
DBETTA*  57.35799      DEGBEES 


HEEL    ANGLE-  11.95908       DEGBEES 

BODDER    AHGLE=  1.  76 328 DEGBEES 
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KTT-  7.50000      KHOTS 

KVA-  10.36289      KHOTS 

SDCCESSPUL   CONVEBGEHCE 
ITER         -         3 

BOAT   SPBED=  5.96515      KHOTS 

LEE8AT    AHGLE=  0.50853       DEGBEES 


DGAHMA=  80.00000      DEGBEES 

DEETTA=  45.45833      DEGBEES 


HEEL    ANGLE*  19.44243      DEGBEES 

BODDEB    fc)IGLE=  3.  1043SDEGBEES 


KVT»  7.50000      KBOTS 

K7A=»  11.17890      KHOTS 

SOCCESSPOL   C0HVE8GENCE 
ITEB         »         4 

BOAT   SPEED=  5.34839      KHOTS 

LEEBAI    ANGLES  1.17645      DEGBEES 


DGAHBA=  60.00000      DEGBBES 

DBETTA=  35.52245      DEGBEES 


BEEL    AHGLE=  24.68011       DEGBEES 

BODDEB    AHGLE=  3. 94630DEGBEES 


K?T=  7.50000      KHOTS 

KVA*  11.17072      KHOTS 

SUCCESSFUL   CONVEBGENCE 
ITEB         *        4 

BOAT   SPEEDS  4.52826      KHOTS 

LEEHAT    ANGLE=  2.44491       DEGBEES 


DGAHHA=  45.00000      DEGBEES 

DEETTA=  28.34309      DEGBEES 


HEEL    AHGLE*  25.90133      DEGREES 

BODDEB    AHGLEs  3. 91684DEGBEES 


K?T=  7.50000      KHOTS 

KVA=  10.53564      KHOTS 

SUCCESSFUL    CONTEBGEHCE 


DGAHHAs  35.00000      DEGBEES 

DEETTA=  24.09880      DEGBEES 
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KVT»  7.50000      KBOTS 

KVA»  10.36289      KNOTS 

SUCCESSFUL    CONVERGENCE 

ITER         -         3 

BOAT   SPBED=  5.96S1S      KNOTS 

LEESAT     ANGLE  =  0.SO853       DEGBEES 


DGAHHA=  80.00000      DEGREES 

DEETTA=  45.45833      DEGBEES 


HEEL    ABGLE=  19.44243      DEGBEES 

BUDDEB    AHGLE=  3. 10435DEGBEES 


KVT*  7.50000      MOTS 

KVA=  11.17890      KHOTS 

SUCCESSFUL   COBVEBGENCE 
ITFB         »         4 

BOAT  SPBED=  5.34839      KNOTS 

LEENAI    ANGLE=  1.17645       DEGBEES 


DGABHA=  60.00000      DEGBEES 

DBETTA=  35.52245      DEGBEES 


HEEL    ABGLE=  24.68011       DEGBEES 

BUDDEB    ABGLE=  3. 94630DEGBEES 


DGAHBA=  45.00000      DEGBEES 

DEETTA=  28.34309      DEGBEES 


KVT=  7.50000      KHOTS 

KV1*  11.17072      KBOTS 

SUCCESSFUL    CONVERGENCE 

ITFB         -        4 

BOAT   SPEED"  4.52826      KBOTS  HEEL    ABGLE*  25.90138      DEGREES 

LEEBAT    AHGLE=  2.44491       DEGBEES  BUDDEB    ANGLE=  3. 91684DEGREES 


K¥T=  7.50000      KBOTS 

KVA»  10.53564       KHOTS 

SUCCESSFUL    CONVERGENCE 


DGABHA=  35.00000      DEGBEES 

DEETTA=  24.09880      DEGBEES 
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ITEH        ■        4 

BOAT   SPEZOs 
LEEBAT    AHGL2= 


3.47308      KHOTS 

4.47  674      DEGREES 


HEEL    AHGL2=  23.64365      DEGBEES 

RODDER    AHGLE=  4. 04634DEGBBES 


KVT=  10.00000       KHOTS 

K7A=  5.16953      HOTS 

SUCCESSFUL   COHVBHGEHCE 
ITER         -         4 

BOAT   SPEED*  4.83  03  7      KHOTS 

LEEBAT    AHGLE*  -0.00401      DEGBEES 


DGAHBA*         180.00000      DEGREES 
DBETTA*         179.99990      DEGREES 


HEEL    AHGLE=  0.00542       DEGBEES 

BUDDEB    AHGLE=  0.  01 243 DEGREES 


KTT=  10.00000       KHOTS 

RVA=  5.48727       KHOTS 

SUCCESSFUL   C0H7EBGEHCE 
ITEH         »         3 

BOAT   SPEED*  5.10574      KHOTS 

LEEBAT    AHGLB*  0.01586      DEGBEES 


DGAHBA=         160.00000      DEGBEES 
DBETTA*         141.44290      DEGBEES 


HEEL    AHGLE*  0.08175      DEGREES 

BUDDEB    AHGLE*  0. 02775DEGHEES 


KVT*  10.00000       KHOTS 

KVA*  6.67725       KHOTS 

SUCCESSFUL   C0H7ERGEHCB 
ITER         =         3 

BOAT   SPEED*  5.85245      KHOTS 

LEEBAT    ANGLE*  0.10  124      DEGBEES 


DGABHA*         140.00000      DEGREES 
DBETTA*         105.70840      DEGREES 


HEEL    AHGLE*  1.40412       DEGREES 

RUDDEB    AHGLE*  0. 17259 DEGBEES 
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KVT=  10.00000       KHOTS 

KVi=  8.79442       KNOTS 

SUCCESSFUL    CONVERGEICE 

ITER*       »         8 

BOAT   SPEED=  6.52654      KHOTS 

LEEJ1T    AHGLE=  0.10446      DEGREES 


DGAHHA=         120.00000      DEGREES 
DBETTA=  79.97887      DEGBEES 


BFEL    AHGLE=  6.55590       DEGBEES 

BUDDER    AHGLE=  1. 10168DEGBEES 


K7T=  10.00000      KHOTS 

KVX=  10.96851       KHOTS 

SUCCESSFUL  COHVEBGEHCE 
ITER         «         ft 

BOAT    SPEED=  6.56510      KHOTS 

LEEBAT    AHGLE=  -0.16  120      DEGBEES 


DGAHHA=  100.00000       DEGBEES 

DBETTA=  63.87666      DEGBEES 


BEEL    ANGLE=  16.25755      DEGBEES 

BUDDEB    AHGLE=  3. 55508DEGBEES 


KVT*  10.00000       KHOTS 

KTA=  12.85525      KHOTS 

SOCCESSFOL  COHVEBGEHCE 
ITER         «         3 

BOAT   SPEED»  6.52192      KHOTS 

LEEBAT    AHGLE=  -0.33663      DEGREES 


DGAHHA=  80.00000      DEGREES 

DBETTA=  50.00259      DEGREES 


BEEL   AHGLE=  27.24298      DEGBEES 

BODDEB    AHGLE=  5. 98752DEGBEES 


KVT=  10.00000       KHOTS 

KVA*  13.80846       KHOTS 

SUCCESSFUL    CONVERGENCE 


DGA?inA=  60.00000       DEGBEES 

DEETTA=  3  8.84161      DEGBEES 
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ITEB         »        5 

BOAT   SPEED* 
LEEBAT    AHGLE= 


5.75414      KNOTS 

0.378U6       DEGBEES 


HEEL    ABGLE=  35.06798      DEGREES 

BODDEB    ANGLE*  7. 60845DEGBEES 


KVT*  10.00000       HOTS 

KVA=  13.68682      KNOTS 

SDCCESSPUL    CONVERGENCE 
ITEB         ■         ft 

BOAT   SPEED*  4.64757      KBOTS 

LEEBAY    ANGLE=  2.15695      DEGBEES 


DGAfl«A=  45.00000      DEGBEES 

DBETTA*  31.10666      DEGBEES 


BEEL    ANGLE*  37.00015       DEGREES 

BDDDEB    ANGLE*  8.  19907DEGBEES 


KVT=  10.00000       KNOTS 

KVA=  12.49910      KBOTS 

SUCCESSFUL    CONVERGENCE 
ITEB         »        7 

BOAT  SPEED*  2.9137  9      KNOTS 

LBEBAI    ANGLE*  7.24246       DEGBEES 


DCAflHA*  35.00000      DEGBEES 

DBETTA*  27.31577      DEGBEES 


BEEL    ANGLE*  32.23439      DEGBEES 

BODDEB    ANGLE*  13. 07657DEGBEBS 
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ITB1Z      ■ 

10 

iphutt  • 

0 

>           - 

4 

EPS1         - 

1.0D-10 

EPS2      *    , 

1.0E-02 

Z0L0(1) ...X0L0(    «) 

1.000000E   01  0.000000E-01  O.00OOO0E-01  0. 000000E-01 

SOLUTIONS    ABE   COBPOTID    POR   THE    FOLLOBING    THOE    BIND    SPEEDS    (KNOTS) 

5.00000         7.50000      10.00000 

SOL0TI0NS    ABE    COMPUTED    FOB    THE    FOLLOBING    TROE    HIND    ANGLES    POB    EACH    BIND    SPEED 

180.00000    160.00000   140.00000    120.00000    100.00000      80.00000      60.00000      45.00000       35.00000 
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L»L=  30.00000 

DISPLACEBEHT    TO    LENGTH'         225.00000 

PRISHATIC=  0.53000 

L/H"  14.76796 

LCB=  53.00000 

KIHEBATIC    VISCOSITY*         0. 1279100B-04 

DENSITY    OF    BATBB*  1.99000 

BEAN    KEEL   C3  0BD=  6.90000 

\ 
BEAN    RODDER    CHORD*  1.87537  s 

BOOGHNESS    AND    CORRELATION    FACTOR*  0.00000 

FAIRBODT    BETTED    SURFACE=         215.25700 

KEEL   WETTED   SURFACE-  63.01086 

BUDDER    BETTED    SURFACE*  19.99725 

TOTAL    BETTED    S0RFACE=         298.26480 

BBL  /LBL=  0.30  276 

SAIL    AREA*         900.00000 

HEELING    HCrtBHT    ABB=  23.99541 

PITCHING    BOBENT   ABB*  26.35417 

DISTANCE    FROH    FP    TO    BODDER    STOCK*  28.49998 

CLOSE   HAD  LED    LIFT    C0EF=  1.25000 

BIG   FORB    DRAG    COEF*  0.40000 

CROSS   SECTIONAL   ABEA   OF    BIG*  26.97595 

SAIL    FBICTIONAL   COEF*  0.04000 

INDUCED    DRAG    COEF   /   LI  PT  C0EF**2*  0.13030 

LATERAL    AREA    OF    HDLL=  121.29310 

TRANSVERSE  AREA    OF    HULL*  37.06689 

DISTANCE    BETBEEN    CE    SAIL    AND    FP*  13.93835 

OFF    BIND    DRAG   COEF*  0.00000 
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OFP    HIHD    LIFT   COEP=  0.00000 

DISTANCE    BETBEEN    CE  HOLl   AMD  PP»  13.93835 

DENSITT   OP  AIB=  0.00238 

HORIZONTAL   DISTANCE  PROA    HAIIflOK    DBAPT   OP    PBOPILE    HEASOBEBENTS 

0.00000         3.00000        6.00000         9.00000      16.50000 

PBOPILE    HEASOBBSENTS 
6.38000        6.25240         1.92985         1.42200         0.00000 

DBAPT    OP    CANOE   BODT    /  TOTAL    DBAPT=  0.3184  1 

HTDRODYNAHIC    LIPT    COEF  =  0.11209 

LAHDA    X    LIPT   /   DBAG=  2.04978 

DISTANCE    FROB    LHL    TO   VEBTICAL   HIDBO    CE=  2.124S4 

***    BODDER    SECTION    *** 

DISTANCE    PBOB   LWL    TO    BOOT  CHOBD=  0.40628 

ROOT   CH0BD=  2.25000 

TIP   CHCRD=  1.50075 

RODDER    SPAN=  5.07766 

SHEEP    ANGLES  0.00000 

ASPECT    BATIO=  5.41509 

BEAN   CHOBD=  1.87537 

DCL  /  DALPHA=  0.07057 

CROSS    FLOW   DRAG    C0EP=  0.56024 

DCB   /  DCL=  0.21654 

POBfl    DRAG   COEF*  0.00650 

BODDER    ABEA=  9.52251 

***    SIGHTING    BOBENT   SECTION    *** 

BETACENTRIC    HEIGHT   /   LHL=  0.11534 

PBOUDE    ND8BEB    RIGHTING    HOHENT   COEP-         0. 60895 1 5s-0« 
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■OK    LIHEAB    BTDBOST1TIC    COEF  = 
DISTANCE    BELOW    LBL    TO    CG= 


0. 3767873 E-03 
0.75000 


***    SAIL    SECTION    **• 

I»  a7. 17540  J=  15.72513 

P=  42.88670  E*  14.29557 

AB   FOB    TRI=  3.00000  AB    nAIM= 

AB   SAIL    PLAN=  2.47280 

CE   SAIL    PLAN    APT    HAST=  0.60041 

PORHABD    OVEBHABG=  2.38720 

PBEEBOABD=  3.71000 


3.00000 
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KVT=  5.00000       KNOTS 

KVA»  2.36043       KNOTS 

SUCCESSFUL   C0H7EBGEHCE 
ITER         ■         5 

BOAT   SPEED=  2.63557      KHOTS 

LEEWAY     ANGLE=  -0.00116       DEGREES 


DGABflA=         180.00000      DEGREES 
DBETTA«         179.99980      DEGRBES 


HEEL    AHGLE=  0.00057       DEGREES 

RODDER    ANGLE=  0. 00 433 DEGREES 


KVT=*  5.00000       KHOTS 

KVA=  2.U913»       KNOTS 

SUCCESSFUL   CONVEHGESCE 
ITER         «         2 

BOAT   SPEED=  2.88579      KHOTS 

LEEBAT    ANGLE?  0.02C97       DEGREES 


DGABBA=  160.00000       DEGREES 

DBETTA=         136.65290      DEGREES 


HEEL   AHGLE=  0.01554      DEGREES 

RUDDER    ANGLE=  0. 00262DEGBEES 


NO   CONVERGENCE 

KTT=  5.00000       KNOTS 

KVA=  3.21880       KHOTS 

SUCCESSFUL   COHTERGEICE 
ITER         »      11 

BOAT   SPEED=  3.65326      KHOTS 

LEEBAT    ANGLE=  0.13603       DEGREES 


DGAHHA=         140.00000      DEGREES 
DBETTA=  93.15175      DEGREES 


HEEL    ANGLE*  0.24055      DEGREES 

RUDDER    ANGLE=  0. 02736 DEGREES 


KVT=  5.00000      KNOTS 

KVA=  4.83768       KNOTS 

SUCCESSFUL   CONVERGENCE 


DGAHBA=         120.00000      DEGREES 
DBETT.\=  63.51920      DEGREES 
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ITEB         »         4 

BOAT   SPEED* 

LEEWAY     ANGLE* 


1.65673      KNOTS 

0.54518      DEGREES 


HEEL    ANGLE=  1.80063       DEGBEES 

BODDEB    ABGLE=  0. 062U0 DEGREES 


KVT*  5.00000       KHOTS 

K?A=  6.53116       KHOTS 

S0CCESSP0L   COHVEBGENCE 
ITEB         *         4 

BOAT   SPEED=  5.15885      KHOTS 

LEESAT    »NGLE=  0.98065      DEGREES 


DGAIJHA=         100.00000      DEGREES 
DBETTA=  48.93198      DEGBEES 


HEEL    ANGLE=  4.21832       DEGBEES 

BODDEB    AMGLE=  0. 20861 DEGBEES 


KVT»  5.00000       KHOTS 

KVA=  7.84539       KHOTS 

SOCCESSPUL   C0H7EBGEHCE 
ITEB         ■         3 

BOAT   SPBED=  5.23940      KHOTS 

LEEUAT     ANGLE=  1.49549       DEGBEES 


DGABMA=  80.00000      DEGBEES 

DBETTA=  38.87585      DEGBEES 


HEEL    ANGLE=  6.78553       DEGBEES 

BODDEB    AHGLE=  0. 37304DEGBEES 


KVT*  5.00000       KHOTS 

K7A=  8.57841       KHOTS 

SOCCESSPOL   COHVERGEHCE 
ITEB         *         3 

BOAT   SPEED=  4.90367      KHOTS 

LEEWAY    AHGLE=  2.18603      DEGBEES 


DGABBA=  60.00000      DEGBEES 

DBETTA=  30.31612      DEGBEES 


HEEL    ANGLE=  8.62737      DEGBEES 

BODDEB    ANGLE=  0. 45298DEGBEES 
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KVT=  5.00000       KNOTS 

KVA=  8.39471       KNOTS 

SUCCESSFUL   CONVERGENCE 
ITEB         ■         4 

BOAT    SPE2D=  4.07831      KNOTS 

LEEH4I    AHGLE=  3.16383       DEGREES 


DGAflHA=  45.00000      OEGBEES 

D£ETTA=  24.90797      OEGBEES 


HEEL    ANGLE-  8.49245      DEGB EES 

BUDDEB    ABGLE=  0. 43054 DBGBEES 


KVT*  5.00000       KNOTS 

KVA=  7.81175      KBOTS 

SUCCESSFUL   CONVEBGENCE 
ITEH         -         4 

BOAT   SPEED=  3.17022      KNOTS 

LEEWAY    AHGLE=  4.64773      DEGREES 


DGAMHA=  35.00000      DEGBEES 

DBETTA=  21.53836      DEGREES 


HEEL    ANGLE=  7.44635       DEGREES 

RUDDER    i»GLE=  0. 38519 DEGB EES 


DGAHHA=         180.00000      DEGREES 
DBETTA=         179.99980      DEGREES 


KVT=  7.50000       KNOTS 

K7A=  3.59100      KNOTS 

SUCCESSFUL   CONVERGENCE 
ITEB        »         3 

BOAT    SPEED=  3.91276      KNOTS  HEEL    AHGLE=  0.00121       DEGREES 

LEEWAY    ANGLE*  0.00031       DEGBEES  RODDER    A»GLE=  0. 00054 DEGB EES 


KVT=  7.50000       KNOTS 

KVA=  3.81778      KNOTS 

SUCCESSFUL   CONVERGENCE 


DGAHHA=  160.00000      DEGBEES 

DBETTA=         137.78660      DEGBEES 
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ITER        ■        2 

BOAT   SPEED" 

LEZSAT    ANGLE= 


4.21828      KHOTS 

0.02009       DEGREES 


HEEL    AHGLE=  0.03173      OEGBBES 

RUDDER    ANGLE*  0. 00429DEGREES 


KVT=  7.50000      KIOTS 

K7A=  4.86734       KIOTS 

SUCCESS  POL   CONVERGENCE 
ITEB        ■         3 

BOAT   SPEED=  5.07021      KIOTS 

LEEWAT    AHGLE=  0.14829       DEGREES 


DGAHBA=  140.00000       DEGREES 

DBETTA«=  97.92107      DEGREES 


BEBL   AHGLE=  0.49961      DEGREES 

RUDDER    ANGLE=  0. 016 19DEGBEES 


KVT=  7.50000      KHOTS 

KVA=  6.84121       KHOTS 

SOCCESSPOL   CONVERGENCE 
ITER        »         3 

BOAT   SPEED=  5.90355      KIOTS 

LEEWAT    A8GLE=  0.51083       DEGREES 


DGAHHA=         120.00000      DEGREES 
DBETTA=  71.69914      DEGREES 


HEEL   ANGLE*  2.91928      DEGREES 

RODDEB    AHGLB=  0. 18288DEGBEES 


K?T=  7.50000      KHOTS 

KVA*  9.05576       KNOTS 

SOCCESSFOL  CONVEBGENCE 
ITEB        «         3 

BOAT   SPE8D=  6.53885      KHOTS 

LEEHAT    ANGLE-  0.95508      DEGREES 


DGABHA=         100.00000      DEGREES 
DBETTA=  54.64867      DEGREES 


HEEL    AHGLE=  7.42339      DEGREES 

BODDEB    ANGLE=  0. 56 130 DEGREES 
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KVT=  7.50000       KJOTS  DGAHfiA=  80.00000       DEGSEES 

KYA=  10.80911      KHOTS  DBETTA=  43.10342      DEGBEES 

SUCCESSFUL    COHVEBGENCE 

ITER         »         3 

BOAT   SPEED=  6.58876      KNOTS  HEEL    ASGLE=  12.24898       DEGBEES 

LEEiAY    ANGLE=  1.51378      DEGBEES  BODDER    AHGLE=  0. 99877 DEGBEES 


K7T*  7.50000      KHOTS 

KVA=  11.74616       KHOTS 

SOCCESSPOL    CONVERGENCE 
ITEB         »         4 

BOAT   SPEED=  6.03699      KHOTS 

LEEWAY    ANGLE=  2.41843       DEGBEES 


DGAHNA=  60.00000      DEGBEES 

DBETTA=  33.57051      DEGBEES 


HEEL    ANGLE=  15.68644       DEGREES 

BODDEB    AHGLE=  1.  26  121  DEGREES 


KVT=  7.50000       KHOTS 

KVA-  11.81018       KHOTS 

SOCCESSPOL   CONVERGEHCB 
ITEB         «         3 

BOAT   SPEED=  5.24840      KHOTS 

LEEHAI    ANGLE'  3.57327      DEGBEES 


DGAHHA=  45.00000      DEGBEES 

DBETTA=  26.68242      DEGBEES 


HEEL    AHGLE*  16.53915      DEGBEES 

BODOER    AHGLE=  1. 20 193DEGBEES 


KVT=  7.50000       KHOTS 

KVA=  11.24568      KHOTS 

SOCCESSPOL   COHVEBGENCE 


DGABHA=  35.00000      DEGBEES 

DBETTA=  22.49055      DEGBEES 
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ITER 


BOAT   SPEED*         '    4.24629      KNOTS 
LEEHAT    ANGLE-  5. 22368      DEGREES 


HEEL    ANGLE*  15.29806      DEGREES 

RODDER    ANGLE*  1. 09039DEGBEES 


KVT=  10.00000       KNOTS 

KVA*  4.96443       KNOTS 

SOCCESSPOL    CONVERGENCE 
ITER         ■         3 

BOAT   SPEED=  5.03554      KNOTS 

LEEHAT    ANGLE=  -0.00046      DEGREES 


DGAHHA*         180.00000      DEGREES 
DBETTA-         179.99980      DEGREES 


HEEL    ANGLE=  0.00228       DEGREES 

RODDER     ANGLES  0.  00276DEGRBES 


K?T=  10.00000       KNOTS 

KVA=  5.32662      KNOTS 

SUCCESSFUL   CONVERGENCE 
ITER         »        3 

BOAT   SPEED*  5.31337      KNOTS 

LEEWAY    ANGLE*  0.01862      DEGREES 


DGAHHA*         160.00000      DEGREES 
DBETTA=         140.05190      DEGREES 


HEEL    ANGLE*  0.04492       DEGREES 

RODDER    ANGLE*  0.  00731  DEGREES 


KTT=  10.00000       KNOTS 

KVA*  6.61337      KNOTS 

SOCCESSPOL    CONVERGENCE 
ITER         *         3 

BOAT   SPEED*  6.10511      KNOTS 

LEEHAT    ANGLE*  0.14843      DEGREES 


DGARBA*         140.00000      DEGREES 
DBETTA*         103.60230      DEGREES 


HEEL    ANGLE*  0.78191       DEGREES 

RODDER    ANGLE*  0. 04948DEGREES 
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KVT»  10.00000      KHOTS 

EVi=  8.88492       KHOTS 

SUCCESSFUL   CONVERGENCE 
ITER         »         9 

BOAT   SPEED=  6.98967      KHOTS 

LEEHAT    1HGL2=  0.45382       DEGREES 


DGAHBA=  120.00000       DEGREES 

DBETTA=  77.08757      DEGREES 


HEEL   AHGLE=  4.09941       DEGREES 

RODDER    AHGLE=  0. 35205DEGHEES 


KTT=  10.00000      KNOTS 

KVA=  11.29239       KHOTS 

SUCCESSFUL  CONVERGEHCE 
ITER         *        8 

BOAT   SPEED=  7.26006      KNOTS 

LEEWAY    ANGLE=  0.89417       DEGREES 


DGAHMA=         100.00000      DEGREES 
DBETTA=  60.70335      DEGREES 


HEEL    AHGLE=  10.27948       DEGREES 

RODDER    AHGLE=  1.  12386 DEGREES 


KVT>  10.00000       KHOTS 

K7A=  13.33439       KHOTS 

SUCCESSFUL   CONVERGENCE 
ITBR         »         3 

BOAT   SPEED=  7.25288      KHOTS 

LEEHAT    AHGLE=  1.48  165      DEGREES 


DGAHBA=  80.00000      DEGREES 

DBETTA=  47.60773      DEGREES 


HEEL    AHGLE=  17.39751       DEGREES 

RODDER    A8GLE=  2.  01724DEGREBS 


KVT=  10.00000      KHOTS 

KVA=  14.65839       KHOTS 

SUCCESSFUL   COHVERGEHCE 


DGAHHA=  60.00000      DEGREES 

DBETTA=  36.21414      DEGREES 
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ITEH        -        5 

BOAT    SPEED* 

LEEWAY    AHGLE= 


6.82657      KHOTS 

2.54576       DEGREES 


HEEL    ANGLE=  23.44096      DEGBEES 

BODDEB    AHGLE=  2. 50608DEGBEES 


KVT*  10.00000       KHOTS 

KTA*  14.70452       KHOTS 

SOCCESSPOL   COHVEBGEHCE 
ITEB         *         4 

BOAT   SPEED*  5.82167       KHOTS 

LEEWAY    AHGLE=  4.21530       DEGBEES 


DGABHA=  45.00000       DEGBEES 

DBETTA=  28.74271      DEGBEES 


HEEL    ANGLE*  24.95070       DEGBEES 

RODDER    AHGLE=  2. 51 346DEGBEES 


KVT*  10.00000       KHOTS 

KYA*  14.17157       KHOTS 

SOCCESSPOL   CONVERGENCE 
ITEB         *         4 

BOAT    SPEED*  4.76666      KNOTS 

LEEUAT    ANGLE=  6.43  894       DEGBEES 


DGABHA=  35.00000      DEGBEES 

DBETTA*  23.87466      DEGBEES 


HEEL    ANGLE*  23.88383       DEGBEES 

RUDDER    ANGLE=  2. 34001  DEGREES 
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Thesis 

M925   Munger 

Prediction  of  sMlino 
boat  performance. 
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